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ABSTRACT 


Ground motion characteristics of the Himalayan earthquakes have 
been studied from the 198 horizontal con^x>nent records suid 99 vertical 
component records from six earthquakes in the region that are now 
available. The characteristics studied are the peak parameters, duration 
of strong shaking, shape of elastic response spectrum and ductility 
reduction factor. It is found that mean A/V ratio is quite high in the 
region ajxi does not vary much with the distance from the shock-source. A 
study of shape of elastic acceleration response spectra revealed that: 
(1) shape varies somewhat insignificantly with the distance, (2) shape 
of spiectra in the intermediate period range highly depends on earthquake 
magnitude, (3) there is a urgent need to revise the shape of response 
spectra recommended in 18:1893-1984, and (4) shape of "standard 
spectra” suggested by AERB for rock sites is compwirable with the 
meeum-plus-one-standard-deviation shape from the data. Both, the ratios 
of peak jiarameters and elaistic response spectra indicate that with 
distance the attenuation of waves of different frequencies is not so 
. different as elsevdiere. Based on the data, shape of ductility reduction 
factor spectrum to be used in design is also proposed. 
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CHAPTERI 

INTRODUCTION 


1.1 GENERAL 

Earthquake-resistant design of structures requires estimate of the 
expected ground motions at the site in future earthquakes. In order to 
develc^ ability to reasonably assess the expected ground motion due to 
future earthquakes, it is essential to study the characteristics of 
strong ground motions recorded in the past earthquakes. 

Design of engineering structures is usually governed by building 
codes and standards of design practices. The reconmendations in these 
codecs are often based on the analysis of records from major earthquakes. 
For example, the seismic design specttixn, specified in the codes, takes 
into consideration characteristics of recorded ground motions. As more 
records become available, these recommendations are updated and improved 
for better eairthquake- resistant design. Unfortunately, until recent 
years, time history records of strong ground motion in the Indian 
earthquakes were almost non-existing. As a result, for instance, the 
shape of design spectrum in the Indian seismic code (IS: 1893 - 1984) is 
based mainly on somie earlier earthquakes in California. 

Since 1986, a number of strong motion records have become available 
from earthquakes in Himalayas. These need to be studied and analysed to 
understard the peculiarities of ground motion in the Himalayan region. 


1.2 SCOPE OF THE STUDT 

In the present work, strong motim instnxnent data from Himalayas, 
recorded during the period 1986-1991, are examined. In all, 198 
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horizontal CGni>anent and 99 vertical conponent records from six 
different earthquakes with nagnitude varying frcm 5.2 to 7.2 and 
^icentral distance in the range of 5 km to 320 km are now available. 
These have been analysed to understaixi the characteristics of 
earthquake ground motion in the Himalayas. The characteristics studied 
are the peak parameters, duration of strong shaking, shape of response 
spectrum, and ductility reduction factor. Results obtained are conpared 
with those from published stvidies based on earthquakes elset^re and 
with currently adopted design practices. 


1.3 ORQANIZATIOM 

Apart from this introductory chapter, the thesis comprises of five 
more chapters to cover the present work. Chapter II gives details about 
the strong motion records used in this study. Chapter III presents the 
characteristics of the data in terms of ratios of peak ground motions 
and duration of strong shaking. Chapter IV presents the study on shape 
of response spectra obtained from the data. Chapter V presents the study 
on ductility reduction factor for elasto-plastic model. Sutimary and 
conclusions from the present work are given in Chapter VI. Appendix A 
gives the details of the records used in the study with regard to 
name of stations triggered in each earthquake, epicentral distances, 
peak ground acceleration, peak ground velocity, peak ground 
displacement, duration of strong shaking (Trifunac and Brady 1975), and 
spectrum intensity (Housner 1959) for 5% damping. 



CHAPTER II 

STRONG NOTION RECORDS 


2.1 DITRODOCriON 

Ground motion during an earthquake is measured by a seismograph 
»diich is essentially a single degree of freedom system with its mass, 
stiffness, and darrping characteristics appropriately designed. 
Instnments with large natural period (flexible oscillator) usually 
record ground displacement; these being flexible are very sensitive and 
can record ev» the slightest motions in ground. Such instruments, 
normally termed as seismographs, are useful for recording very small 
earthquakes or for recording large earthquakes taking place at large 
distance away from the iiistrument. However, these go off-scale in case 
of a strong shaking, which is of prinary interest for engineering 
applications. On the other hand, instnments with a stiff oscillator 
(large natural frequency) measure ground acceleration and are not 
sensitive to weak ground motions. Such instnments, often termed as 
strong motion accelerographs, provide the realistic levels of ground 
acceleration that an engineering structure may be subjected to during a 
strong shaking. In normal operations, these instnments lie in a passive 
state, and get triggered, to operate, when the ground motion exceeds a 
threshold value set for the instrument. 

Hie acceleration versus time records obtained from such instnments 
are called accelerograms. Usually, an accelerogram is in the form of an 
analog record, such as a photographic film. However, in recent years, 
digital instnments have also become available. The analog records 
obtained from the accelerograph films have to be digitized at closely 
spaced time intervals. The digital time history, obtained directly or 



4 


from digitizing the analog record, is then corrected through appropriate 
filtering to control the digitization errors and baseline distortion; 
also certain accelerograph transducer corrections are applied. The 
records thus processed are termed as the "corrected” accelerograms. The 
corrected accelerogram can then be integrated with respect to time to 
cA)tain a velocity verst;is time record, which can be fxarther integrated to 
dt)tain a displacement versus time record. 

2.2 CHARACmiSATION OF STRONG GROUND MOTION 

Pig. 2.1 shows some typical accelerograms obtained from epicentral 
region of different earthquakes. A visxial comparison of these records 
suggests that the characteristics of earthquake ground motion may vary 
significantly with regard to the peak acceleration, nunnber of 
significant peaks in the time record, and duration of shaking. 

Fig. 2. 2 shows the acceleration, velocity, and displacement time 
histories derived from a typical accelerogram. Not surprisingly, the 
integration of record changes its frequency content; vMle the 
acceleration record has large high frequency content, the displacement 
record shows motion of high period, and the velocity record shows an 
intermediate frequency content. A structure is expected to achieve 
resonance (or multifold increase in response) when frequency of ground 
motion matches with natural frequency of the structure. Hence, 
structures with different natural period are expected to respond 
differently in the same earthquake. 

The earthquake ground motion, vdrLch is somewhat random with time, 
can be characterized by several parameters d^ending on the application. 
Hie parameters of most interest are the peak parameters (peak ground 
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acceleraticai, velocity, or displacement), duration of strong shaking, 
and its frequency content. Hie frequency content is somewhat more 
difficult to quantify and is usually represented through a Fourier 
spectrum or a response spectnm. 

Traditionally, in the eaurthquake- resistant design, peak ground 
acceleration has been the single most frequently used parameter for 
describing the severity of ground motion. A rigid structxire subjected to 
a given eairthquake ground motion experiences a total earthquake force 
equal to its mass times the peak ground acceleration. Thus, peak ground 
acceleration gives a good physical feel of the ground motion and is a 
convenient rule of thunnb for scaling different accelerograms. The 
difficulty with peak ground acceleration is that, sometimes, it gives a 
very misleading indication of the earthquake ground motion. For 
instance, an earthq[uake of rather low magnitude (say M = 4 or 5) may 
give a rather large peak ground acceleration, even though the overall 
motion may not be very damaging to structures. 

The peak ground velocity is associated with waves of intermediate 
frequencies (mostly between 0.3 to 3.0 Hz), which is also the range in 
which most civil engineering structures lie. Haice, peak ground velocity 
is a much better indicator of the damage potential of a ground motion 
than the peak ground acceleration. Peak ground displacement is 
associated with waves of low frequencies (0.1 to 0.3 Hz); in this 
frequency range very few structures lie, hence this parameter is not 
very useful for engineering applications. However, peak grotmd 
displacement is an inportant parameter for seismological studies, 
because, amplitudes of ground displacement are propagated more 
coherently than amplitudes of ground velocity and acceleration. This due 
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to the fact that their low-frequency spectral composition is not very 
saisitive to scattering by small geologic inhomogeneities, ilso, peak 
ground displacement gives a somendiat better indication of the fault 
displacement. 

2.3 STRONG MOTION INSm»lENTATION IN HIMALASAS 

Ihe arc-shaped Himalayan belt extending from Pamirs in erstwhile 
USSR to Arkans in Burma has long been regarded as a highly seismic area 
in the world. The Hinalayan belt is believed to have formed due to 
collision of Indo-Australian and Eurasian plates in past several ndllian 
years. The area has experienced several destructive earthquakes of 
magnitude greater than 8.0. The notable among these earthquakes are the 
Assam earthquake of 1897 (M = 8.7), the Kangra earthquake of 1905 (M = 
8.6), the Bihar-Nepal earthquake of 1934 (M = 8.4), and the Assana-Tibet 
earthquake of 1950 (M = 8.7) (Big 2.3). Nuorterous earthquakes of 

nagnit\xies ranging between 6.0 and 7.0 have also occurred in the region. 

An international workshop on strong motion earthquake instrument 
arrays held in Honolulu, Hawaii, in May 1978 (Iwan 1978) Identified the 
area around Shillong in the north-eastern India as the highest priority 
site for deployment of a strong motion array. Subsequently, three strong 
motion arrays, listed below, have been deployed in the Hinalayan region: 

(a) An array of fifty strong motion accelerographs (SMA-1 of nv/s 

Kinemetrics, USA) is located in north-west Himalayas in the state of 

Himachal Pradesh. This array is identified as Kangra array. 

(b) An array of fortyfive SMA-1 is located in north-east Himalayas in 
the states of Assam and Meghalaya. This array is identified as 
Shillong array. 
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(c) An array ccnprising of forty SMA-1 instnments is located in the 
central Himalayas in the state of Uttar Pradesh. This is identified 
as Uttar Pradesh array. 

Fig. 2.4 shoMS the location of these arrays. There are also plans 
to instrument the Himalayas of Bihar-Nepal border region. Thus, the 
Himalayan zone is now better instrumented to record the strong-motion 
data from future earthquakes. 

Besides the above three arrays, the region also has a nunnber of 
indigenoijs strong motion accelerographs and structural response 
recorders (BQ 88-16 1988, BQ 93-16 1993); these have also provided 
records of strong ground motion. However, the data from these 
instirunents has not been considered in the present stxidy. 

2.4 SraONG MOTICW DATA USED IN THE PRESENT STODT 

The Kangra array registered an event of magnitude 5.5 on April 26, 
1986. This is the only significant event recorded so far by this array. 
The Siillong array has registered foxn: events during 1986-1988, while 
the Uttar Pradesh array has registered one event in October 1991. The 
details of these events with regard to their date, magnitude, maximum 
peak ground acceleration, and the range of epi central distance are 
presented in Table 2.1. This table also gives the abbreviaticxi used in 
rest of this study for describing different events (e.g., "Sept” for the 
earthquake of September 10, 1986). The magnitudes indicated in table are 
as per estimates of the United States Geological Survey (USGS); these 
differ somewhat from those by the India Meteorological Department (IMD). 
Thus, the present study is based on 99 sets of strong motion records 
(198 horizontal components and 99 vertical conponents) obtained from six 
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earthquakes (magnitude ranging from 5.2 to 7.2) that have ocxwrred in 
three different regions of the Himalayas. The epicentral distance for 
these records varies from about 5 km to 320 km. Detailed description on 
site geology of stations is not available. However, ChandrasheMiaran and 
Das (1990) indicate that the stations are located on "reasonably firm 
ground", "soil", or on "soft sedimentary rock". In the absence of more 
specific information about the local soil conditions, all stations are 
assianed to have almost similar ground conditions, i.e., the records have 
not been classified on the basis of soil conditions. Moreover, idiile 
conparing with other data, the stations are assuned to be at rock sites. 

Figures 2.5 to 2.10 show the location of instruments triggered 
during the different earthquake events. Also shown in these figures is 
the location of epicenter. Pigs. 2.5 to 2.9 show three locations of 
epicentre: El is the estimated locaticxi by the IMD, E2 is that by the 
USC^, and E3 is the location suggested by Chandrashekharan and Das 
(1990) based on records from the array. For "April" (Fig. 2.5), revised 
estimate of epicentre by IMD, E4, is also shown. For each of the events. 
Appendix A gives the name of the stations triggered and three different 
estimates of epicentral distance for each station. Also presented in 
Appendix is the peak ground acceleration, peak ground velocity, peak 
ground displacement, spectrum intensity (Housner 1959) for 5% damping 
and duration of strong shaking (Trifunac and Brady 1975) for each of the 
records used in the study. The different estimates of epicontre in these 
earthquakes vary by as much as 68 km for "Sept" and 85 km for "Feb". 
Such large variaticxi rules-out the possibility of using epicentral 
distance as a criteria for classifying the records for low, medium, and 
large epicentral distances. Moreover, for most of these events, neither 
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an indication of rupt\nre zone (based on aftershock data) nor the 
isoseismals are available. Thus, the distance of recording station from 
the zone of energy release also cannot be iised as a criteria for 
classifying the records. However, in this study, epicentral distance 
refers to the epicentral location suggested by Chanshdchran and Das 
(1990) based on the analysis of strong motion records. The location of 
^i centre for the "Oct” event (Fig. 2.10) is fairly accurate. 

To accoimt for the variation in epicentral distance, each set of 
records has been classified into one of the three groups on the basis of 
the valtie of larger horizontal peak ground acceleration (PGR). Records 
with Isurger horizontal peak ground acceleration less than 0.05g are 
classified into the low PGA group, vhereas those with greater than O.lg 
PGA are classified into the high PGA group; records with larger 
horizontal peak ground acceleration between these two limits are 
classified into the mediizn PGA group. Thxas, once the leurger horizontal 
PGA from a set falls in one groxjp, all the three ccmponants of that set 
are included in that groxip. Hie limits of PGA for classification were 
chosen so that in each grov:®) a reasonable number of accelerograms could 
be included. Table 2.2 shows the siirroary of nunber of records included 
in each of the groups from different events. 



Table 2.1 Summary of the earthquake events considered in this study. 


Rrray 

region 

Date of 
event 

Abbreviated 

as 

M 

NS 

R (km) 

POR (g) 

Min. 

Max. 

Hbrz 

Vert 

Kangra 

Rpril 26, 
1986 

April 

5.5 

9 

5 

29 

0.250 

0.082 

Shillong 

Septenber 10, 
1986 

Sept 

5.2 

12 

15 

181 

0.138 

0.062 

MsmKMgm 

May 


14 

46 

179 

0-086 

0.005 

February 6, 
1988 

Pdt> 

5.8 

18 

35 

176 

0.114 

0.103 

August 6, 

1988 

Aug 

7.2 

33 

100 

324 

0.343 

0.180 

Uttar 

Pradesh 

October 20, 
1991 

Oct 

6.6 

13 

25 

151 

0.310 

0.294 


M : Magnitude of earthquake 

NS : Nunber of stations triggered 

R : l^icentral distance (approximate) 

PGR : Peak ground acceleration 

Horz : Horizontal connpcxient 

Vert : Vertical ccniKxvent 


Table 2.2 Sixmary of nunnber of records in the three PGR groups. 


Array 

region 

Event 

High PGR 
group 

Medium PGR 
group 

Low PGR 
group 

Horz 

Vert 

Horz 

Vert 

Horz 

Vert 



10 

5 

4 

2 

4 

2 



6 

3 

8 

4 

10 

5 


■H 

0 

12 

6 

16 

8 

CUXLl XOilfJf 

Feb 

2 

1 

12 

6 

22 

11 

BSHU 

28 

14 

32 

16 

6 

3 

Uttar 

Pradesh 

Oct 

8 

6 

12 

4 

6 

3 
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Fig. 2.2 Accel eraticm, velocity, and displaconent time 
histories for the SOOE corponent of El Centro 
earthuake of 1940 (Mohraz and Elghadamsi 1989). 
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CHAPTER III 


RATIOS (F PEAK C®OUND MOTION AND 
DURATION OF STRONG SHAKING 


3.1 INTRODOCTIOII 

Given an estimate of horizontal peak ground acceleration (PGA) at a 
given site, the values of peak ground velocity (PGV) and displacement 
(PGD) are often obtained by using average values of ratio of peak 
parameters. Two ratios, A/V and AD/V^ (A: peak ground acceleration, V: 
peak grotind velocity, D: peak ground displacement) are most often used 
for this purpose. Similarly, vertical peak ground acceleration is 
usxially estimated by using an average value of the ratio •' 

vertical peak ground acceleration, : horizontal peak ground 

acceleration). A nixiber of studies are available vMch indicate the 
average value of these ratios for a given type of site condition, based 
on groxjnd motions recorded elsewhere in the world (e.g. , Mohraz et al. 
1972, Mohraz 1976, Mohraz 1991). 

Seismological studies show that attenuation of seismic waves is 
frequency dependent; high frequency waves attenuate more rapidly with 
distance than moderate or low frequency waves. Since, peak ground 
acceleration, peak ground velocity, and peak groimd displacement are 
associated with waves of different frequencies, these ratios of peak 
parameters happen to be good indicators of ground motion characteristics 
at different sites. For example, ground motions e3q)erienced near an 
earthquake source are expected to have higher value of VV than ground 
motions at a large distance from the source of seismic energy release. 
This is so because the high frequency waves (and hence the PGA) 
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attenuate faster with disteuice than the intermediate frequency waves 
(and hence the PGV) . 

It can be seen that for hanimic oscillations, AD/V = 1.0. With 
distance, the ground motion tends to become harmonic with 
ever-increasing period; hence as epicentral distance increases, the 
AD/V ratio must reduce and approach to xmity (Newmark and Rosenblueth 
1971). Also, certain response spectrum characteristics can be correlated 
to the AD/V^ ratio. 

This chapter presents a study of the peak parameter ratios and 
duration of strong shaking from strong motion data of the six Himalayan 
earthqueikes . 

3.2 REVIEW OP TOE PAST STODIES 
3.2.1 A/V Ratio 

Mohraz et al. (1972) conducted the first study on V/A ratio. Tliey 
obtained the average values of V/A from a data set, mostly consisting of 
records obtained from the San Fernando earthquake of 1971. They 
recotrmended the value of V/A ratio as 48 in/sec/g (A/V ratio as 8.1 
sec ^) for firm ground, and 32 in/sec/g to 36 in/sec/g (A/V ratio in the 
range 11.7 sec ^ to 12.1 sec ^) for rock sites. 

Mohraz (1976) studied the variation of V/A ratio with different 
site conditions using 54 earthquake records obtained in 16 seismic 
events in California. The records were classified into four categories 
depending on the geological conditions at site. For each site category, 
the ratios were obtained for three groups; the first group consisted of 
horizontal component records with larger PGA, the second of horizontal 
component records with smaller PGA, and the third of vertical component 
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records. Table 3.1 sunreirises these results. It can be seen that fl/V 
ratio is higher for rock sites than for allvtvium. Also, K/V ratio for 
horizontal caTi>onent groi® with larger PGA is higher than that for the 
lower PGA group. A/V values for the vertical coni>anent are close to 
those for horizontal component groi;p with larger PGA. 

Variation of V/A ratio (and hence A/V ratio) with distance from the 
nflptured surface (R) has been studied by Mchraz (1991) from records of 
the Lana Prieta earthquake of 1989 (Table 3.2). It is seen that at rock 
sites VV ratio for near-field recoxrds (R < 20 km) can be more than 
twice that for mid-field records (50 km > R > 20 km). At allxxviixn sites, 
the difference is about 30 %. However, as the site distance increases 
beyond 50 km, A/V ratio shows a slight increase. 

Using 36 records from rock sites, Zhu et al . (1988) studied the 
effect of ground motions with low, normal, and high A/V ratio on the 
damage of single degree of freedom stiffness-degrading systems. They 
demonstrated that if a system is s\ji>jected to ground motions with same 
PGA as that of the design spectrum, the ground motions with low A/V 
ratio can cause much more significant peak inelastic deformatioi, 
stiffness deterioration, and hysteretic energy dissipation than is the 
case for ground motions with high A/V ratio. This is consistent with the 
general expectation that given two ground motions with the same PGA, the 
motion with higher PGV will be more damaging to the structures. 

Tso et al. (r992) have examined the engineering implications of A/V 
ratio on the ground motim characteristics. A data set consisting of 45 
records obtained from rock sites during 23 events having nagnitude 
ranging between 5.2 to 8.1 was used. The records were classified as low, 
nomal, or high A/V when this ratio was less than 0.8 g/m/s ( = 7.9 
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sec~^), between 0.8 g/m/s and 1.2 g/m/s ( = 11.8 sec and more than 

1.2 g/ra/s. The stvdy concluded the following: 

- Ground motions in the vicinity of small or moderate earthquake usually 
have high A/V ratio, whereas those distant from large earthquakes have 
low A/V (Fig. 3.1). 

- Ground motions with high A/V ratio usually have seismic energy in the 
high frequency range, whereas those with low A/V ratio have aoergy in 
low frequency range. It was seen that ordinates in acceleration 
spectra for periods longer than about 0.5 sec for low A/V records 
could be more than twice as high as those for high A/V group (Fig. 
3.2). 

- Peak ground velocity correlates well with the spectrum intensity 
(Houfflier 1959) over earthquake records having drastically different 
A/V ratio. 

3.2.2 AD/V^ Ratio 

2 

The AD/V ratio was first introduced by Esteva (1969). The ratio is 

primarily used to estimate PGD from known PGA and PGV. He suggested an 

2 

enpirical equation for estimating AD/V for a known epicentral distance 
(R in km): 

AD . T . 400 

rO.6 

Newmark and Rosenblueth (1971) suggest that choice of AD/V to evaluate 
PGD is preferable as compared to D/A or D/V ratios because the 
integrated velocities and displacements are sensitive to baseline 
adjustment and large errors may accumulate dxxring the integration. By 
selecting AD/V , the error which may arise is distributed among the 
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ground motion parameters. According to Newmark and Rosenblueth, for most 

2 

earthquakes of practical interest 2VD/V ranges ' from approximately 5 to 
15. 

Mohraz et al. (1972) recottinend an "adequately conseirvative'' value 

of 6.0 for AD/V^ ratio. Further, they observed that AD/V^ is a measure 

of sharpness or flatness of pseudo velocity spectrum, i.e., hig^ber the 
2 

value of AD/V ratio, flatter will be the pse\jdo velocity spectrum. 

2 

Table 3.1 clearly shows that AD/V is higher for rock sites than for the 
alluvium; however, this is not so obvioxis from the Lena Prieta 
earthquake data (Table 3.2). Also, Table 3.2 does not show any specific 
trend in the value of AD/V ratio with site distance from the ruptured 
surface. 

3.2.3 Ay/^^ Ratio 

Table 3.3 shows the mean and 84.1 percentile values of the vertical 
to larger horizontal PGA ratio (A^/Ag) eJotained by Mohraz (1976). For 
rock sites the mean value is around 0.52 while the 84.1 percentile value 
is around 0.69. He observes that 'The ratio of the vertical to larger 
peak horizontal acceleration of 2/3, which is currently enployed in 
seismic design corresponds to 84.1 percentile ratio and not the mean. 
Although this value is conservative, its use has been justified to 
account for variations greater than the median and mean". The Indian 
code (18:1893-1984) recommends that the ratio of vertical to horizontal 
seismic coefficient be taken as 0.5. 

3.3 RESULTS ERCM THE DATA 

For each of the six events, mean and mean-plus-standard deviation 
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2 

values of the h/V and AD/V ratios for high, medinn, and low PGR grov^js 
have been connputed; these are sumarized in Table 3.4. The ratio for 
vertical PGA to the larger horizontal PGA (A^/Ag) for differaat events 
is given in Table 3.4. The trends for the Himalayan earthquakes from 
these tables are discussed here. 

3.3.1 A/V Ratio 

Table 3.4 clearly shows that the VV ratio for the Indian 
earthquakes is unusually hi^. For the horizontal components, the mean 
value for different events ranges from 13.3 sec ^ to 32.5 sec ^ with an 
overall mean value of 20.2 sec These values are to be compared with 
mean values of 12 sec ^ to 15 sec ^ for rock sites suggested by 
different studies of Mohraz based cm data for elsewhere. For instance, 
Mohraz (1972) suggests the value as around 12 sec Table 3.1 indicates 
a value of 14.5 sec and the data of the Loma Prieta earthquake (Table 
3.2) indicates a mean value of 14.8 sec ^ for near-field (and 
significantly lower values for site distance exceeding 20 km). 
Similarly, the data studied by Tso et al. (1992) pertains to ^icaitral 
distance ranging from 4 km to 379 km; of all these the data pertaining 
to "high” A/V ratio of 11.7 sec or more comes from motions recorded 
within 27 km of the epicentre. 

Table 3.4 also shows that A/V ratio is generally higher for the 
high PGA group data than for the meditm or low PGA group data; i.e., the 
expected trend of decrease in A/V ratio with increase in the epicentral 
distance is generally seen. The only exception to this are the 
horizontal ground motion records of event "Oct” and the vertical grovmd 
motion records of event "April”. However, the variation in A/V ratio 
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with distance is not so significant for the Indian data as it is for the 
other eaurthquakes . For instance, in the present case, the epi central 
distance ranges from about 5 km to 320 km; for such distances, A/V ratio 
for earthquakes elsewhere could vary by more than 100%, which is not the 
case in the Himalayan data. 

Generally, the A/V ratio for vertical ccrponent of ground motion is 
cotrparable to that for the horizontal components. 

On the whole, a value of around 20 sec”^ could be considered as an 
appropriate mean value and 12 sec~^ as the 
mean-minus-one-standard-deviation value for the Himalayan earthquake 
records. This is much higher than the values recommended in the 
literature based on records elsetdiere. This has very serious 
implications for earthquake-resistant design of important structures in 
the Himalayas. This means that the earthquake ground motion in the 
Indian Himalayas are associated with (a) unusually low value of ground 
velocity and hence are far less damaging, (b) unusually hig^i value of 
ground acceleration and hence are more damaging (particularly to the 
stiffer structures), or (c) a combinaticai of the possibilities (a) and 
(b) above. To resolve this issue, systematic studies on attenuation 
relationships for the Himalayan earthquakes are required for which far 
more authentic information on the earthquake source parameters is 
needed. In the meanwhile, for seismic risk studies in India, a great 
deal of cauticn is required while using the empirical relationships 
(relating PGA and/or PGV with the earthquake magnitude and ^icentral 
distance) developed for data from elsewhere. 

Table 3.5 shows the mean value of A/V ratio obtained in each of the 
six Himalayan earthquakes. The results are tabulated in increasing order 
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of earthquake magnitude. It is readily apparent from the table that mean 
A/V valtoe decreases with the increase in magnitude. This is probably 
because, for a given PGR, the value PGV increases with the increase in 
earthquake magnitude. For the three earthquake magnitude ranges; (1) 
small magnitude (M < 5.5), (2) moderate magnitude (5.5 s M s 6.5), and 
(3) large magnitude (M > 6.5), the mean A/V value are also presorted in 
the same table. 

3.3.2 AD/V^ Ratio 

Table 3.4 also shows the valxies of mean AD/V^ for the three groups 
in each event. AD/V^ for vertical component is higher than it is for the 
horizontal conponent; this is caisistent with the data analysed by 
Mohraz (1976) and seen in Table 3.1. On the whole, the mean value of 
AD/V ratio for the horizontal component is 4.6; this is lower than the 
value of around 7.0 for rock sites in Table 3.1, but higher than 2.3 to 
3.8 based on the Loma Prieta data (Table 3.2). The data for "i^ril” has 
one record with unusually high value of AD/V (equal to 85.6) which 
could be due to a typographical error in the paper by Chandrasekaran and 
Das (1990). This gives vtnusually high values of the mean and 
mean-plus-one-standard-deviation for low PGR group of this evait. Haace, 
in arriving at the overall mean aiad mean-plus-one-standard-deviation 
values for the vertical components (5.2 and 8.8, respectively), this 
single value has been disregarded. The data of Table 3.4 show no clear 
trend in the variation of AD/V with distance, except for the event 
"Aug". 


3.3.3 fiy/Ajj Ratio 
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A statistical svtimary of the ratio of in the three PGA 
groups is presented in Table 3.4. iMs table shows that, except the 
event ”Oct”, the mean value of ratio is around 0.5 and 
mean-plus-one-standard-deviation value is arovmd 0.70. This event shows 
a rather strong vertical conoponent of ground motion with a mean value of 
A^lig as 0.66 and mean-plxjs-one-standard-deviation value of 0.9. For the 
entire data, the mean value of Ay/A^ ' is 0.51 and the 
mean-plus-one-standard-deviation value is 0.70. These values are similar 
to those reported by Mohraz (1976) for data from elsewhere (Table 3.3). 
Thus, to be on the conservative side, the design for vertical 
acceleration should be based on ratio as around 0.67 or 0.75 as 
against the value of 0.5 presently being used in the Indian code. 


3.4 DORATION OF STRONG SHAKING 
3.4.1 Definition 

The total duration of a complete accelerogram represents the 
duration during which a earthquake gromd acceleration exceeded a 
threshold or critical acceleration value set for the instrumait. 
However, this duration may not represent the duration dxiring which the 
structure was subjected to severe shaking. In the past, several 
procedures have been suggested for extracting a strong motion segment 
from an accelerogram; each of these give different value of the duration 
as seen in Fig. 3.3. As yet, no universally accepted definition for 
duration of strong motion is available. 

Bolt (1969) and Page et al . (1972) have suggested a 'bracketed 

dtiration"; this is the duration between the first and the last 
acceleration peaks exceeding a specified value (usually 0.05g). Since 
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PGA decreases with increase in disteuice from epicenter, tMs definition 
leads to a decrease in the duration of strong motion with an increase in 
distance from the source. 

Trifunac and Brady (1975) define the duration of strong shaking as 
the time duration during which the middle 90 percent (i.e., 5% to 95%) 
contribution of the "acceleration intensity" (defined as jA^dt) takes 
place. They observe that the average duration for soil is approximately 
10-12 sec longer than that for rock, and that it increases at the rate 
of 1.0 - 1.5 sec for every 10 km increase in distance from the 
shock-source. 

McCann and Shah (1979) define the duration of strong ground motion 
based on the concept that the strong motion part can be defined as the 
period during which the RMS, or average rate of energy, is highest 
relative to the rest of the record. In this method, first the root mean 
sguare (RMS) value io ovaluntod at e.ich of tJir> digit ir-rnl in tho 
accelerogram, this is referred to as cumulative root mean square 
function (CRF). Rate of change of CRF during total duration of the 
record is then examined to give two cut-off points. The final cut-off 
time is obtained at a time where rate of change of CRF becomes negative 
and remains so for remainder of the record. The initial cut-off time is 
obtained in the same manner, except that the search is performed in the 
reverse direction of time scale starting from the end of the record and 
continued till beginning of the record. 

3.4.2 Coni>arisaQ of Different Definitions 

A graphical coirparison of duration as pet the Trifunac and Brady 
(1975) and McCann and Shall (1979) definitions is shown in Fig. 3.4 for 
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sane records fron the Californian earthquakes. This figxire shows that 
duration by the McCann and Shah method is shorter than that by the 
definition of Trifunac and Brady (1975). On the other hand. Table 3.7 
shows the average duration of records from the Lana Prieta earthqixake 
based on three different definitions. This table indicates that the 
definition by McCann and Shah gives duration longer than that by the 
Trifunac and Brady method. This table also shows that the duration is 
generally higher for alluvixin sites than for the rock sites. 

3.4.3 Variation' in Duration irith Distanca and with A/V Ratio 

Tso et al. (1992) demonstrate that the records with high JV/V ratio 
have shorter durations of strong shaking than those with intermediate 
A/V ratio, whereas the records with low A/V ratio have longer durations. 
They explain this in terms of epicentral distance as, "Larger the 
magnitude of an earthquake, the longer the duration of strong ground 
motion if the distance from the epicerxtre remains const 2 mit. If the 
magnitude of the earthquake is kept constant, duration increases with an 
increase in epicentral distance. Since the records with high h/V ratios 
were generated in the proximity of small or moderate earthquakes, 
whereas those with low A/V ratios were obtained at large distances from 
large or moderate earthquakes, the correlation of the A/V ratio with the 
duration of strong shaking can be expected." 


3.5 DURATION OF STRONG SHAKING FOR. THE DATA 

The duration of strong shaking has been calculated for each of 
thel98 horizontal records based on the Trifunac and Brady (1975) 
definition (Appendix A). The relationship between the duration and A/V 
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ratio for each of the events is shown in Fig. 3.5. The points for 
records belonging to different PGA groups (low, medium, and high) are 
shown with a different notation. The plots of Fig. 3.5 show no definite 
relationship or trend between duration and A/V or between duration and 
PGA groxap (i.e., epi central distance), though such relationships have 
been reported in literature for studies based on data from other 
earthquakes (Tso et al . 1992) . 

This can also be seen in Table 3.8 which shows the mean duration of 
strong shaking in the three PGA groxj^s for the six Himalayan 
earthquakes. For "kpril" and ”Fd3", the trend of increasing duration 
with decrease in PGA (i.e., increase in distance from shock-source) is 
quite clear. But, exactly opposite trend is observed in the event "Aug”, 
where high PGA group records have greatest duratio of strong shaking. 
For events "Sept" and ”Oct”, the Table shows the largest duration of 
shaking in the medium PGA group records. However, Fig. 3.5 clearly 
confirms the fact that larger the earthquake magnitude, grater is 
duration of of strong shaking. For *'sept” (M = 5.2) and "April" (M = 
5.5), duration of most of records is less than 10 sec, while for events 
"Oct" (M = 6.6) and "Aug" ( M = 7.2), duration of most of the records is 
more 10 sec. 
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Table 3.1 Surtnary of V/A and AD/V^ ratios for earthquakes in 
California (Mohraz 1976) . 


Site category 

Group 

Mean 

V/A 

(in/s)/g 

Approx, 
mean A/V 

(Sec"^) 

AD/V^ 

Mean 

84.1 

percentile 

Rock 

m 

27.0 

14.5 

6.9 

11.0 

s 

30,0 

13.0 

IHQQI 

11.2 

B 

31.0 

12.6 

7.6 

11.8 

Less than 30 ft of 
allixviunn underlain 
by rock 

B 

37.0 

10.6 

5.2 

7.7 

s 

43.0 

mm 

5.2 

8.2 

B 

37.0 

10.6 

8.5 

13.3 

30 - 200 ft of 
alltivixxn underlain 
by rock 

B 

33.0 

11.8 

5.6 

7.8 

s 

41.0 

9.5 

4.3 

mm 

B 

33.0 

11.8 

9.1 

13.7 

Alluviun 

m 

51.0 

BH 

4.3 

6.0 

s 

62.0 

6.3 

B 

4.4 

Bl 

51.0 

IQH 

5. A 

7.0 


L : Horizontal corponents with the larger peak ground acceleration 
S : Horizontal components with the smaller peak groijnd acceleration 
V : Vertical components 
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Table 3.2 Sumtnary of mean V/A and AD/V ratios for different 
source to site distances in Lana Prieta earthquake 
of 1989 (Mohraz 1991). 


Category 

Rock 

Alluvium 

V/A 

(in/s)/g 

Approx. 

A/V 

(sec ^) 

AD/V^ 

V/A 

(in/s)/g 

Approx . 
A/V 

(sec ) 

AD/V^ 

Near-field 
(R < 20 km) 

26.0 

15.0 

3.8 

40.0 

9.8 

2.9 

Mid-field 

(50 km > R > 20 km) 

63.0 

6.2 

2.3 

53.0 


2.9 


54.0 

7.2 

2.8 

48.0 

8.1 

2.0 


R : Distance of recording station fron ruptured surface 


Table 3.3 Sunmary of the ratio of peak vertical to larger 
peak horizontal acceleration (Mohraz 1976) . 


Site category 

Mean 

84.1 

Percentile 

Rock 

0.52 

0.69 

Less than 30 ft of alluvium underlain 
by rock 

0.49 

0.62 

30 - 200 ft of alluvium underlain 
by rock 

0.46 

0.66 

Alluvium 

0.45 

0.61 
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1^1 e 3.4 Sumary of A/V and AD/V^ ratios for the six aurthquake 
events in Himalayas. 


Event 

PGA 

Group 

Horizontal 

coni>onents 

PGA 

Groifl? 

Vertical 



AD/V^ 


AD/V^ 

m 

m - o 

m 

m + o 

m 

m - Oj 

m 

m + <7 


mm 

21.1 

14.4 

2.9 

4.2 

H (5) 

23.2 

17.9 ’ 

5.6 

.7.8 

April 

M (4) 

21.4 

15.7 

4.2 

B 

M (2) 

14.3 

13.8 

3.0 

3.5 


L (4) 

15.7 

7.2 

4.2 

IQQI 

L (2) 

22.7 

HIQQII 

M.8 

89.5 


H (6) 

32.5 

25.8 

3.5 

B 

H (3) 

36.0 

30.1 

9.2 

13.6 

Sq?t- 

M (8) 

31.0 

20.0 

5.9 

10.8 

M (4) 

26.8 

22.8 

10.1 

12.5 


L(10) 

26.5 

20.1 

B 

8.6 

L (5) 

22.9 

13.6 

6.9 

9.4 


H (0) 








B 


May 

M(12) 

23.7 

18.3 

6.5 

9.9 

M (6) 

19.5 

11.8 

f.3 

11.2 


L(16) 

17.3 

12.6 

5.0 

6.5 

I. (8) 

15.2 

12.5 

4.8 

6.5 


H (2) 

22.9 

21.3 

5.0 

B 

H (1) 

28.7 

M 

9.2 


Feb 

M(12) 

22.6 

17.4 

HQQI 

B 

M (6) 

25.1 

17.7 

8.1 

11.6 


1.(22) 

18.7 

12.7 

5.2 

8.1 

L(ll) 

15.2 

9.2 

4.6 

6.2 


H(28) 

23.3 

17.1 


5.9 

H(14) 

25.8 

18.5 

6.8 

9.4 

Aug 

M(32) 

16.7 

9.8 

1 3.6 

5.6 

M(16) 

16.3 

7.8 

4.4 

7.5 


L (6) 

13.3 

B; 

2.6 

1 

3.5 

L (3) 

12.8 

8.9 

3.2 



H (8) 

14.3 

12.2 

IQQ 


H (6) 

15.0 

10.3 

2.9 

3.6 

Oct 

M(12) 

19.1 

13.0 

B 

6.8 

M (4) 

14.4 

■n: 

2.9 

4.3 


L (6) 

14.7 

11.3 

3.2 

4.6 

L (3) 

10.7 

9.0' 

2.9 

3.2 

All 

(198) 

20.2 

12.5 

4.6 

8.0 

(99) 

19.6 

1 

10.8 

5.2 

8.8 


* : Number of records 

H : High PGA group (PGA > 0.1 g) 

M : Medium PGA group (0.1 > PGA > 0.05 g) 

L : Low PGA group (PGA < 0.05 g) 

m : Mean 

m + <7 : Mean-plus-one-standard-deviation 
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Table 3.6 Suttmary of peak vertical to larger horizontal 
acceleration ratio for the data used in the 
present study. 


Event 

PGA 

Group 

V®H 

V*H ° 

f each event 

- - 


m 

m + a 

m 

m + a 


asi 

0.34 

0.46 



April 

M (2) 

0.57 

0.76 

0.50 

0.74 


L (2) 

0.81 

0.97 




H (3) 

0.45 

0.47 



Sept 

M (4) 

0.41 

0.52 

0.44 

0.51 


L (5) 

0.40 

0.60 




H (0) 



.. 



May 

M (6) 

0.39 

0.52 

0.51 

0.72 . 


L (8) 

O.W 

0.81 




H (1) 

0.89 

0.89 



Feb 

M (6) 

0.41 

0.47 

0.50 

0.72 


h(ll) 

0.52 

0.77 




H(14) 

0.44 

0.55 



Aug 

M(16) 

0.49 

0.63 

0.47 

0.61 


L (3) 

0.59 

0.67 

1 



H (4) 

0.85 

1.05 



Oct 

M (6) 

m 

0.77 

0.65 

0.90 


L (3) 

0.53 

0.80 



All 

(99) 



0.51 

0.70 


* : Munnber of records 

H : High PGA group (PGA > 0.1 g) 

M : Medium PGA group (0.1 > PGA > 0.05 g) 

L ; Low PGA group (PGA < 0.05 g) 

m : Mean 

m + o ; Mean-plus-’one-standard-deviation 

















































































38 


Table 3.7 Average duration of strong motion (in sec) for accelerograms 
recorded in the Loma Prieta earthquake of 1989 (Mohraz 1991). 
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Table 3.8 Sunnary of duration of strcaig 
shaking for the records obtained 
from the six Himalayan earthquakes. 


Event Magnitude 
of event 




Duration 


mean 


3.48 


M (4) 


L (4) I 10.55 


H (6) 


M (8) 


L(10) 


H (0) 


8.94 


mean 


5.60 



8.33 


M(12) 16.27 16.28 

L(16) 17.07 


H (2) 10.51 

M(12) 13.11 16.48 

L(22) 18.86 


H(28) 26.10 

M(32) 22.83 23.65 

L (6) 16.62 


H (8) 11.69 

M(12) 17.00 14.47 

L (6) 13.10 


: Nunber of records 
: Hi^ PGA group (PGA > 0.1 g) 

: Mediun PGA group (0.1 > PGA > 0.05 g) 
: Low PGA group (PGA < 0.05 g) 
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EPICENTRAL DISTANCE (KM) 

Fig. 3.1 Relationship between earthquake magnitude (M) and 

epi central distance (R) for low, intermediate, and 
high A/V records (Tso et al. 1992). 



PERIOD (SECOND) 


Mean 5% damped elcistic acceleration response spectra 
for the three A/V groups and whole ensemble of 
records scaled to a ccninon PGA of Ig 
(Tso et al. 1992). 


Fig. 3.2 





Accelerafton, tn/sec/sec 
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Time, sec. 

Fig- 3.3 Ccrrparison of strong motion duration for the S69E 

conponent of the Taft, earthquake of July 21, 1982 
using different procedures 
(Mc^az and Elghadamsi 1989). 



Duration by 

Trifvmac and Brady (1975) 


Fig- 3.4 Corrparison of duration of strong shaking from 
Trifunac-Brady (1975) and McCann-Shah (1979) 
definitions (McCann and Shah 1979). 




April (M = 5.5) 


Sept (M = 5.2) 


0 .30 ,40 50 60 
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May (M=5.7) 


0 50 60 


Feb (M = 5.8) 


0 50 60 


0 50 60 


Aug (M=7.2) 


Oct (M = 6.6) 


4 

A ° A tp6 D ° 

a 


o o'" , 

, 0 A I 


A A 

A 


O (S> ^ 


6 1 0 20 30 _,40 50 60 

aA (s®c ; 


10 20 30 ,40 50 60 
A/V (sec’’) 

DDDDD High PGA group 
AAA A A Medium PGA group 
ooooo Low PGA group 


Fig. Relationship between duration and A/V ratio for the 

six Himalayan events. 












CHAPTER IV 

ELASTIC RESPONSE SPECTRUM 


4.1 nUROTOCTIOIf 

Response spectrum is a very usefxd way to characterize strong 
ground motion from the view point of engineering applications. Respite 
spectnm of a ground motion is a plot showing maximum response induced 
by that ground motion in single degree of freedom (SDOF) systems of 
different natural frequency (or nattiral period), but for a given value 
of damping. The response here can be any quantity of interest siich as 
the displacement, velocity, or acceleration. For engineering purposes, 
the following response quantities are most conmcxily plotted (Fig. 4.1): 

(a) Maximum relative displacement (displacement spectnxn) 

SD(«,C) = l«(t)|^ (4.1) 

t 

(b) Maximum relative velocity (velocity spectrum) 

SV(«,C) = l^iCt)!^ (4.2) 

t 

(c) Mzudnum absolute acceleration (acceleration spectrum) 

sa(<.>,C) = |ii(t) + %(t)l^ (4.3) 

t 

(d) Pseudo spectral velocity (pseudo-velocity spectrum) 

PSV(u,C) = » SD(«,C) ' (4.4) 

(e) Pseudo spectral acceleration (pseudo-acceleration spectmn) 

PSA(«,C) = SD(«,C) 

fdiere u(t), Ti(t), ii(t) ace relative displacement, velocity, and 
acceleration, respectively, at instant t for a SDCF oscillator of 
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natural frequaicy « and danping C; and ^ groimd acceleration. 

A response spectnm represents the coinbined influence of anftlitude 
of ground motion, frequency content, and to some extent the duration of 
ground shaking, on different structures (Seed aixi Idriss 1982). Most 
importantly, a response spectrum provides a most convenient way of 
evaluating the maxiiiun acceleration, and thus naximutn inertia force, 
developed in a structure subjected to the given base motion, 
irrespective of fdiether structure behaves as a SDOP system or MDCF 
system. Also, since a SDOF system responds mostly to the energy in the 
ground motion near its am natural frequency, the response spectnin 
characterizes the energy in ground motion. 

4.1.1 Tripartite Representation of Response l^pectrum 

It is obvious from Eqs. (4.4) and (4.5) that SD, PSV, and PSA are 
directly related and provide essentially the same information; hence the 
corresponcLing response spectra can be presented on a tripartite or a 
four-way logarithnic plot by a single curve for each value of damping. A 
typical response spectnin cm a tripartite-plot is shown in Pig. 4.2, 
whidh is essentially a plot of PSV versus T on log-log scale with two 
additional log scales at 45** to period scale for reading SD and PSA. Ihe 
parameters PSA and PSV have certain characteristics that are of 
practical interest. The PSV is close to SV for high frequencies, and 
practically equal for intermediate frequencies, but different for low 
frequencies. PSA is practically equal to SA for damping levels 
encountered in most aigineering applications. 


4.1.2 Design Spectnin 
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Response spectnxn, being a plot of maxiinin response (of a SDCF 
oscillator with varying natural frequaacy) to a given gromd motion 

f 

record, is a description of that particular ground motion. Hence, its 
use is suitable for analysis of a structure subjected to that ground 
motion; but it may not be suitable for the purposes of design (e.g., 
Anderson 1989). Furthermore, a single earthquake record has a particular 
frequency content fdrich gives rise to the jagged, sawtooth appearance of 
peaks and valleys as shown in Big. 4.2. When such a spectrun is used for 
design, stnictures with slis^tly different natural periods may reflect 
vastly different maxinum response and hence the design parameters. To 
overcome this problem^ the conc^t of smooth-shaped response spectnxn 
has been introduced for design. 

A usual method to obtain a smooth-shaped spectrum is through 
normalization and averaging of several response spectra obtained from 
records with cornnon characteristics. Amngst the suggested normalization 
procedures for response spectnm, the most popular is normalization to 
unit peak ground acceleration; at times the normalisation is done with 
respect to other parameters such as peak ground velocity or 
displacement. Principles of statistics are then used to create a 
smooth-shaped spechnxn. The mean or median shape of response spectra so 
obtained is generally xised for eairthquake-resistwt design of normal 
building structures. Use of such a spectrun implies that for a given 
value of peak ground acceleration, there is about 50% probability that 
the maximun acceleration of the SDOP oscillator will be exceeded. High 
risk stnictures may be deigned to higher level such as the 
mean-plus -one-standard-deviation shape, vhich implies that the 
probability of exceedance is only about 16%. To account for the effect 
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of magnitude and distance of different likely events in the vicinity of 
a major project, site-specific design spectnm is sonetimes obtained by 

I 

drawing an envelqp of more than one spectra (e.g. , Fig. 4.3). 

4.1.3 The Actual Design l^>ectrunn 

The above discussion pertains to the shape of design spectnin for a 
given notmalised acceleration, velocity, or displacement. The actual 
design spectnin is obtained Iqr multiplying the entire spectnin shape 
with the design ground motion parameter, for instance the maximum ground 
acceleration (zero period acceleration, ZPA). The magnitude of this 
design parameter has to be decided considering the expected ground 
motion at the site, the factors of safety aiKl load factors to be used 
in design, the ductility to be provided in the structure, the darrping 
characteristics, and the procedure to be used for estimation of natural 
period for the structiire (e.g., Housner and Jennings 1982). 

4.1.4 The Chjective of the Study 

A statistical stuiy of the response spectnin shape has been carried 
out for the records obtained from the Himalayan earthquakes, with a view 
to see if there are any appreciable differences in the ground motion 
characteristics of the Himalayan earthquakes vis-a-vis earthquakes 
elsewhere. A conparison has also been made with the spectnin shapes 
currently in design use. 

4.2 REVIEH OP PAST STUDIES ON DESIGN SPECTPm 
4.2.1 Housner’s Spectnm 

While Biot (1941) initiated the application of response spectnin 
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to earthquake engineering in 1941/ it was not tntil 1959 that Housner 
develc^>ed the first smooth design spectra (Housner 1959). Iliese (Fig. 
4.4} were obtained tising eight strong motion records frcxn four 
earthquakes recorded in California. The curves were arbitrarily scaled 
to 0.125g at zero period. The Indian code (18:1893-1984) has adopted 
the shape given by Housner with minor modifications (Fig 4.5). 

4.2.2 Nemark and Hall ^pectrun 

Another set of design spectra, used extensively in the niaclear 
power plant design in late sixties and early seventies, were suggested 
by Newmark and Hall (1969). They proposed a procedure for ifstimating 
site-independent design spectnxn, which was based on the fact that 
response spectrum over certain frequency ranges is related by an 
amplification factor to the peak value of ground acceleration, velocity 
and displacement. To appreciate this, consider the response spectrum of 
Fig. 4.2 and compare it with the peak ground displacement, peak ground 
velocity, and peak ground acceleration values marked on the same plot. 
It is seen that at small frequencies, the maxiimm relative displacement 
is large, whereas pseudo-acceleratiai is very small. It large 
frequencies, the relative displacement is extremely small, whereas the 
psexjdo-acceleration is relatively large. At intermediate fregoencies, 
the pseudo-velocity is substantially larger than at either end of the 
spectrum. For this reason, three regions are usually identified in a 
response spectrum as low-frequency or displacement region, the 
intermediate-frequency or velocity region, and the high-freqaency or 
acceleration region. In each region, the corresponding ground nation is 
amplified the most. Fig. 4.2 also shows that at small frequencies (0.05 
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Hz or less) the spectral displacement approaches the peak ground 
displacement, indicating that for very flexible systems the maximum 
relative displacement in a simple oscillator is equal to the ground 
displacement. At large frequencies (30 Hz or more), the 

pseudo-acceleration aiqproaches the peak ground acceleration, indicating 
that for rigid systems the absolute acceleration of the mass is same as 
that of the ground. Newmark and Hall suggested that the three 
amplifications (displacement, velocity, and acceleration) are constant 
in their respective regions and can be determined statistically. Tables 
4.1 and 4.2 give the amplification factors sxiggested by NeMonark and Hall 
for horizontal ground motions for median and 84 percentile cases. Thus, 
in this method, first the estimated values of peak ground acceleration, 
velocity, and displacement for the site are plotted on a tripartite 
logari Untie paper; then using the amplification factors, for the desired 
percentage of critical damping, the peak ground motion values are 
amplified to give a smooth design spectnzn for the site. Their spectra 
for ’’standard" earthquake with peak ground acceleration, velocity and 
displacement of l.Og, 121.8 cnv'sec and 91.4 cm, respectively, are shown 
in Fig. 4.6. It can be noted from the figure that curves for different 
damping values approach the peak ground acceleration at different 
period. Newmark and Hall observed that thoiigh the amplification factors 
for horizontal and vertical spectra remain unchanged, the vertical peak 
ground accelerations are about two-thirds of the horizontal peak ground 
accelerations. 

4.2.3 Study by Blxjne and Mohraz 

After San Fernando earthquake of 1971, two independent studies, one 
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by Blxme et al. (1972) aiid the other by Mohraz et al. (1972) , were 
carried out with the objective of developing recornmendations for 
horizontal and vertical design response spectra for nuclear facilities. 
These studies were based on a much larger number of earthquake 
ground-motion records than iJised by the Newmark and Hall study. Ihe 
results of these studies formed the basis for the guidelines contained 
in Nuclear Regulatory Coranreission (NRC) Regulatory Guide 1.6 (1973). The 
technique of representation was very similar to the one proposed by 
Newmark and Hall, except that these stxidies recognized that the 
amplifications in the three regions need not be constant always and may 
vary linearly. To define horizontal response spectra, the peak 
horizontal ground acceleration and displacement levels are first 
established. The peak horizontal ground displacement was considered to 
be proportional to the peak horizontal ground acceleration and was fixed 
at 91.4 cm for an acceleration of Ig (for a standard earthquake). The 
bounds of each of spectnxn are established by five line segments, 
similar to the Newmark and Hall procedure. The amplification factors and 
the response spectra for "standard" earthquake, for horizontal ground 
motion, as adopted by NRC are shown in Fig. 4.7. The vertical response 
spectra were constructed in a manner similar to horizontal response 
spectra, but three differences were incorporated into the procedure: (1) 
the frequency control points were shifted slightly, (2) the 
amplification factors were different, and (3) the value of peak 
horizontal acceleration is used as the initial reference value. The 
amplification factors and the response spectra for vertical ground 
motion are shown in Fig. 4.8. The NRC spectra was also used extensively 
for seismic design of many non-nuclear facilities. 
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4.2.4 £^;>ectnin for Different Soil Conditions 

In 1976, two sttidies, one by Seed et al. (1976) and the other by 
Mohraz (1976), considered the influence of soil conditions on the shape 
of response spectnm. The study by Seed et al. used a total of 104 
horizontal conponents of earthquake records with peak ground 
acceleration greater than 0,05g from 23 different earthquakes. The 
records were divided into four categories based on soil condition at 
recording station. The mean and mean-plus-one-standard-deviatioo spectra 
for the foxxr soil types, obtained by Seed et al. are seen in Pigs. 4.9 
and 4.10, respectively. The figures show that for periods greater than 
about 0.4 sec, the amplifications for rock are substantially lower than 
those for soft to medium clay and for deep cohesionless soil. The AERB 
safety guide (S-11 1990) recommends standard response spectra for rode 
and soil sites based on the work of seed et al. (1976); Pig. 4.11 shows 
the shape recommended by AERB for rock sites. 

4.2.5 Study by Mohraz 

Mohraz (1976) considered a total of 162 records and divided these 
into four soil categories. To exclude the influence of the horizontal 
component of records with smaller peak ground acceleration on 
amplification, he divided the response spectra for each soil category 
into three sets: horizontal components with the larger of the two peak 
ground accelerations, horizontal components with the smaller of the two 
peak ground accelerations, and vertical componofits. Figure 4.12 shows 
the mean acceleration spectra for different soil conditions for the 
horizontal components with larger of the two peak ground accelerations. 
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t is seen that the amplification for allwiutn extends over a larger 
xequency region than the amplifications for the other three soil 
ategories. Also, for short periods, the spectral ordinates for ail^Ivium 
ire lower than for the other soil types, whereas for intermediate and 
Long periods, these are higher. 

Hius, local soil conditions at the site significantly alt«r the 
frequency content in the ground motion and hence the response spectra. 
In the recent y^rs, this has lead to the adoption of different design 
spectra for different soil conditions. Fig. 4.13 shows one such spectra 


adopted by OBC (1991). 

Mohraz (1978) studied the influence of earthquake magnitucje on 
acceleration amplification for alluvium. As seen in Fig. 4.14, 
acceleration amplifications are larger in earthquake of magnitude 
between 6 and 7 than for those with magnitude between 5 and 6. 

Mohraz (1991) examined the influence of source-to-distance on 
spectral shapes for the records obtained from the Lome Prieta earthquake 
of 1989. The accelerograms obtained from rock and alluvium were divided 
into three subgrotqjs based on distance of the recording station from the 
ground projection of the ruptured surface. The groups were labeled as 
"near-field” for records with distance less than 20 km, 'W.d-field" for 
records with distance between 20 and 50 km, and "far-field" for records 
with distance greater than 50 km. Fig. 4.15(a) shows the spectrid shapes 
for the three subgroups obtained for 5% damping using horizontal 
components of accelerograms recorded cn rock sites. The figure shows 
that the amplification for the near-field is substantially smaller than 
the amplification for mid- or far-field for periods greater thm about 
0.5 sec. For periods less than 0.5 sec, however, the amplification for 
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the near-field is greater. Pig 4, 15(b) shows the results of similar 
exercise for records obtained from alluvium. The figure indicates that, 
though distance influences the spectral shapes, the influence is not as 
pronounced as it is for rock. 

4.3 MEmOD OP ANRLYSIS 

Hie data set consisted of 198 horizontal and 99 vertical conponent 
records, obtained from the six events referred to in chapter II. The 
ordinates of the pseudo-acceleration response spectra (PSA) are computed 
by the Newmark-beta numerical method of solution (e.g.. Bathe and Wilscm 
1978). The ordinates are evaluated for a period range of 0.0075 sec to 
3.0 sec at an interval of 0.0075 sec, thus generating about 400 points 
in this range. The computations are carried out for 0, 0.5, 2, 5, 10, 
and 20 percent of critical darnping. Since the peak ground acceleration 
for various accelerograms differ, the computed responses cannot be used 
in statistical study without normalization. In this stixiy, normalization 
of spectra by peak ground acceleration is adopted. Thus, each of the 
pseudo-acceleration spectra is normalized by respective PGA, giving zero 
period ordinates of unity. 

The normalized spectra in each of the events are separated for 
horizontal and vertical components. These are then averaged to obtain 
mean and mean-plus-one-standard-deviation spectral shape for horizontal 
and vertical components representative for individual events. Similarly, 
a mean and mean-plus-one-standard-deviation spectra for horizontal and 
vertical component, considering all the records from 6 events is also 
evaluated. The average spectra are smoothened by adopting a 
five-point-averaging scheme. The "smooth-shapes" thus obtained are then 
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ccnpared with the spectral shapes ciurr^tly in use. 


4.4 RESULTS 

Tlie mean pseudo-acceleration spectra for 6 events for danpinsf 
values of 0, 0.5, 2, 5, 10, and 20 percent of critical danping for 
horizontal and vertical components, are shown in Pig. 4.16. 

4.4.1 General shape 

A visual compariscai of spectra in the Pig. 4.16 suggests that in 
all the events, for a given damping value, there is a build up in 
acceleration amplification frcm ijinity at zero period to a peak value, 
which usually occurs in the period range of 0.1 to 0.3 sec. Por periods 
longer than 0.3 sec the acceleration amplification drops down. However, 
for periods greater than about 2.0 sec the acceleration amplification is 
quite low and remains almost constant. 

Por low damping values (less than 5%), the spectral shape and the 
value of peak amplification in different events vary significantly. 
However, for higher damping values, the peaks of spectra in different 
events are nearly same and the variation in the shape is also 
comparatively less. For instance, in the event "Sept” (M = 5.2) the peak 
anplification for 0% aiKi 5% damping are about 4.8 and 2.4, respectively, 
while in "Aug" (M = 7.2) the corresponding amplifications are about 9.4 
and 2.6, respectively. Thus, a significant variation to the extent of 
200%, is observed in zero damping spectra while the difference is 
insignificant for 5% damping. Though the shape of spectra differs from 
event to event, the peak acceleration amplification in the six events 
for 5, 10, and 20 percent of critical damping is nearly same and is 
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about 2.6, 2.0, and 1.4, respectively . 

4.4.2 Variation within the Data 

Spectra in the "Sept" event drops down sixMenly at period of about 
0.3 sec. Hus could be due to its low magnitude. ’iMay" and "E'eb" are of 
almost same magnitude, their spectral shapes are in general similar, 
except for some variation in maxinun amplification for danping less than 
5% of critical . "Aug" is of the largest magnitude, and in this 
earthquake most stations were at large epicentral distances (greater 
than 150 km) . The spectra for this event show comparatively (with 
respect to spectra for the other earthquakes) more energy in 
intermediate period range (0.4-1. 2 sec); this is expected because as 
distance from the shock-source increases, the high period waves 
dondnate. In case of "April", though the epicentral distances involved 
were quite low (less than 30 km), the ordinates in intermediate period 
range show a substantial amplificaticn. This is unexpected, since 
generally seismic waves in near-field (low epicentral distance) have 
high frequency content. "Oct" spectra are unique in the sense that only 
these show a flat top as compared to sharp tips observed in spectral 
shapes of the other events. 

4.4.3 Variation in Shape of Spectrum with Earthquake Magnitude 
Considering the large variation in amplifications for low values of 

danping, and the fact that in the design of ordinary buildings a danping 
value of 5% is usually adopted, the variation in spectral shape with 
earthquake magnitude has been studied for 5% danping. Fig. 4.17 shows 
the mean spectra for the six events for 5% danping. It is seen that 
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acceleration, an^plifications in the period range of 0.3 to 1.5 sec are 
highly dependent on earthquake magnitude. Larger the magnitude, greater 
is the amplification in this range. For instance, at natural period of 
0.8 sec, the acceleration amplification is 0.3 for the event "Sept” (M = 
5.2) and it is 1.0 for the event "Aug" (M = 7-2), giving a variation of 
as much as 500%. For periods less than 0.2 sec the variation is 
negligible, while it is not very significant for periods longer than 1.5 
sec. This figure also suggests that one could obtain somewhat 
representative shapes of the response spectrum based on earthquake 
magnitude. Fig. 4.18 shows shapes suggested by this data for rock-sites 
for (a) smal 1 -magnitude earthquakes (M < 5.5), (b) moderate-size 
earthquakes ( 5.5 ^ M s 6.5), and (c} large earthquakes (M > 6.5). Ihese 
have been obtained from (a) the mean shape of event "Sept" (M = 5.2), 
(b) average of the mean shapes of events "l^ril" (M = 5.5), '^May" (M - 
5.7), and "Feb" (M = 5.8), and (c) average of the mean shapes of events 
"Oct" (M = 6.6) and "Aug" (M = 7.2). 

The variation of acceleration amplification in the intermediate 
period range with earthquake magnitude can also be anticipated from the 
A/V ratio. As seen in table 3.5, in the low-magnitude earthquakes, the 
value of A/V is high, i.e., for the given value of PGA, the value of PGV 
is low. In the intermediate period range (the velocity region), the PGV 
governs the spectral ordinat^, and hence, low acceleration 
amplification are obtained in this range. 

Not nany stijdies seem to have been done on the shape of 
pseudo-acceleration spectrum with earthquake magnitude. The study by 
Mohraz (1978) (Fig. 4.14) shows this variation but for alluvium sites. 
While the study by Mohraz based on data from elsewhere shows that ground 
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motions of large magnitude earthquakes have significantly larger values 
of acceleration amplification in all the three regions (acceleration, 
velocity, and displacement), the data of the Himalayan earthquake (Fig. 
4.17) suggests that this is so only for the velocity region. 

It is obvious that when choosing spectrum shape for preparing 
site-specific design spectrum, the earthquake magnitude should be 
explicitly considered, particularly in the intermediate period range; 
this is the period range in which a very large nunfiber of industrial 
structures lie. Presently, such a consideration seems to be absent in 
developing the design criteria for important projects in India. 

4.4.4 What does Mean £i>ectra from Data Represent ? 

Mean horizontal pseudo-acceleration spectra, derived by averaging 
all 198 horizontal component spectra is shown in Pig. 4.20. Comparison 
of shape of mean spectra from all horizontal records with spectral 
shapes from "Oct" (M = 6.6) and "Feb" (M = 5.8), shown in Fig. 4.21, 
indicates that mean shape (from all horizontal records) lies in between 
the mean spectra for two evaits, suggesting that mean spectral shape 
from all records is a representative of earthquake of magnitude of about 
6.2 for moderate distances (80-130 km). Incidentally, the average of 
magnitudes of events studied is 6.0. 

4.5 VARIATION IN SPECTRAL SHAPES WITH DISTANCE 

To study the variation in response spectrum shape with epicoatral 
distance, mean response spectra for high, medium and low horizontal PGA 
groups are plotted for different events in Pig. 4.19. The shapes for 
events "Sept", "May", and "Oct" do not show significant variation with 
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distance even though the peaks in the acceleraticm region vary somev^t. 
For "April", the ordinates at large values of natural period are higher 
for the low PGA group, but in this event epi central distance for all the 
stations lie within 30 km. For event "F^", the shape for high PGA group 
is significantly different from that for the other two groups, but in 
this case only two records (i.e., one station) falls in the high PGA 
group and therefore this comparison may not be very general. In the 
intermediate period range, the events "Feb" and "Aug" show significant 
variation for different PGA groups; the high PGA groi:®) (low epi central 
distance) has much lower values of acceleration amplification in this 
period range as is also expected. 

Thus, in general, the Himalayan data does not show a significant 
variation in the spectral shape with epicentral distance; this is in 
contrast with the data from elsewhere. For instance, ccirpare the 
variation for rock sites from the Loma Prieta earthquake (Fig. 4.15) 
with the event "Oct" (Fig. 4.19). While the Loma Prieta earthquake shows 
a very significant variation even within the epicontral distance of 80 
km, the data of "Oct" pertaining to the epicentral distances of 25-150 
km shows not much variation. Only for event "Aug", where large 
epicentral distances (100-320 km) are involved, there is significant 
difference in spectral shape with distance, that too only in the 

velocity region. 

On the whole, it could be concluded that the Himalayan data shows 
somewhat insignificant variation in spectral shape with epicentral 
distance, except in the intermediate period range for the very large 
distances (say more than about 150 km). This is in line with the 
observation made in Chapter III that A/V ratio for the Himalayan 
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earthquake does not vary as significantly with distance as the data from 
elsewhere would suggest. Iherefore, it is obvious that the geological 
featxxres in the Himalayas are such that the attenuation is not 
significantly different for waves of different frequenci^. 

An examination of mean spectral shapes in the three PGA groups 
together with mean A/V of the corresponding groi; 5 j (Table 3.4) reveals 
that spectra with high A/V value have least amplification in moderate 
period range, whereas for spectra in low A/V group just the i^posite is 
true. This again demonstrates that A/V is a good indicator of frequency 
content of ground motion. 

4.6 MEAN-PLUS-ONE-STANDAPD-DEVIATION SPEC?nRtlM 

While the mean spectnm shape is often reccrmended for design of 
buildings, important facilities such as nuclear power plants are usually 
designed on the basis of mean-plus-one-standard-deviation shape. This 
ensures that the variation in spectra ordinates for a givai value of PGA 
is accounted for in an adequately conservative manner. In the case of 
Himalayan data, where many parameters such as local soil conditions, 
epi central distance, and earthquake magnitude are not very reliably 
known, the mean-plus-one-standard-deviation shape has the added 
significance such a shape based on all the data accoiaits for variation 
in these factors. 

Fig. 4.22 shows the mean-plus-one-standard-deviation spectra for 
the six events. In general, shape of mean-plus-one-standard-deviation 
spectra is very similar to shape of mean spectra (Fig. 4.16). Similar to 
mean spectra for the three range of earthquake magnitude (Fig. 4.21), 
nean-plus-one-standard-deviation spectra are also obtained. Fig. 4.23 
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shows the average representative shapes of the spectra for large 
nagnitude (M > 6.5), moderate nagnitude (5.5 ^ M ^ 6.5) and low 
magnitude (M < 5.5) earthquakes. Fig. 4.23 reveals that the 

airplification in intermediate period range is highly dependent on 
earthquake magnitude, fts nagnitude increases, the spectral ordiixates 
also increases and the increase could be as high as 500%. 

Figure 4.24 shows the mean-plus-one-standard-deviation spectra 
evaluated considering all the horizontal component records; these are 
the shapes that may be adopted for major projects in the Himalayas. 

Figure 4.25 shows the coefficient of variation (COV) of response 
spectra for 0, 5, and 20 percent damping, considering all the horizontal 
ccmponent records. The figure suggests that (XN decreases with increase 
in danping. Also, <XN is highly dependent on the natxiral period. The 
value of crv for 5% danping is as low as 0.17 at period of 0.1 sec and 
increases to a value of around 0.8 at the period 1.2 sec. This is 
because the scaling of spectra has been perfonned with respect to the 

PGA. 


4.7 COMPARISON OF SPECTRAL SHAPES 

The design spectrum adopted by 18:1893-1984 is compared with the 
the mean horizontal spectral shape for 5% danping in the six Himalayan 
earthquakes (Fig. 4.26). It is readily apparent from the figure that 
irrespective of magnitude of earthquake, the IS code spectra 
underestinates the acceleration amplification in period range of 0.1 to 
0.35 sec while overestimates in the period range greater than 0.35 sec. 
The mean spectra of event "Aug" (M = 7.2) is closest to IS spectra for 
periods greater than 0.5 sec. In this event, records were obtained from 
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large epi central distances (100-320 km); this suggests that for long 
period range, the IS spectra actually represents large-magnitude, 
far-field spectra. IS code spectrum shape in the low-period range needs 
to be revised upwards in light of this data. As will also be seen in the 
next chapter, this is the period range in whicdi ductility is somewhat 
less effective in reducing the seismic force on a structiire. 

4.7.1 Mean Spectrum with the IS code and UBC (1991) Spectrum 

A coirparisan of the mean spectral shape based on all 198 records 
with the design spectra of 13:1893-1984 and UBC-1991 (rock site) has 
been shown in Fig. 4.27. The shape given in UBC-1991 appears quite 
conservative as compared to the Himalayan data. Note however that 
UBC-1991 does not specify a value of damping for this spectrum shape; it 
is assumed here to be 5%. 

4.7.2 Mean-Plus-One-Standard-Deviaticn Spectra with the AEEQB spectra 

Fig. 4.28 shows comparison of mean-plus-one-standaml-deviation 
spectral shapes from the data with the ’’standard spectra” suggested by 
AERB for rock sites. It is seen that except for zero damping, the AEE® 
spectra match very well with the mean-plus-one-standard-deviation 
spectra obtained from the data. 

4.8 STUOT ON AMPLIFICATION FACTORS 

Representation of spectra in the characteristics trapezoidal form 
as done in the previous works by Newmark and Hall (1969), Mokaz et al. 
(1972), Blume et al.(1972), and Mohraz (1976) has also been attempted in 
the present stxidy. For this the ordinates of PSV, PSA, and SD spectra 
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are evaluated using Newmark-beta method (Bathe and Wilson 1978). The 
ordinates were computed for a set of 200 periods, distributed between 
0.003 sec to about 30.0 sec. The intervals of periods are chosen such 
that they get unifonnly distributed on a log scale. 

From the log-log plots of mean and mean-plus-one-standard-deviation 
normalised PSA, PSV, and SD spectra in each of the events, the most 
aniJlified region is identified. The period limits betweai which these 
amplifications occur are recognised aus amplification region of that 
spectra. For instance, the most amplified region in a PSA spectra is 
recognised as acceleration region. The spectrum in the most amplified 
region is then represented by a straight line, which in most cases is a 
good approximation. This represoitation is illustrated in Figs. 4.29 and 
4.30, which show mean PSA and PSV spectra, respectively, for "Oct" 
event. Also shown in those figures are the straight line representation 
of the spectra in the acceleration and velocity regions, respectively. 

Figure 4.30 also suggests that spectra in the displacanent region 
(which is mostly between 3 to 10 sec) cannot be represented as a 
straight line as has been done in the previous studies. In fact the 
shape of spectra in this region is unexpected. Further investigations 
are necessary for the representation of spectra in this regicm. Hence, 
in the present study amplification factors only for the acceleration and 
velocity regions are evaluated; this is sufficient for most of the 
engineering structures. 

The spectra in acceleration regicm. are obtained as straight lines 
parallel to period axis, that is, constant amplification (constant with 
respect to period) of peak ground acceleration is obtained. However, in 
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inclined straight lines (e.g.. Fig. 4.30), indicating that in this 
period range the amplification of PG7 is linearly varying (on log-log 
plots) and it decreases as period increases. For each of the events 
considered, the mean acceleration and velocity amplification factors are 
presented in Table 4.3 and Table 4.4. For velocity region, the spectra 
being linearly varying, amplification factors at the ends of the 
velocity region are presented. The Tables also indicate the period 
limits of the acceleration and velocity region as seen for these 
events . 

Table 4.3 -reveals that for low damping values (less than 5%) the 
acceleration amplification factor varies for the six events, however 
for higher damping values this factor is quite consistent for the events 
considered. Also, the period limits defining the acceleration region do 
not vary significantly for different events. Generally, acceleration region 
of the Himalayan data is betweei 0.12 sec to 0.26 sec. Table 4.4 
indicates large variation in amplification factors in the period range 
defining the velocity region. Though the lower period end of velocity 
region is somewhat consistent in the six events, the higher period end 
varies significantly. This also gives rise to large variation in 
amplification factor at the higher period end of velocity region. The 
velocity region for the Himalayan data lies between 0.26 to 2.6 sec. 

Acceleration and velocity amplification factors evaluated 
considering all the horizontal component records, are ccnpared with the 
amplification factors reported in previous studies by Newmark and Hall 
(1969), Mohraz (1976) and those recontmended in NRC guide 1.6 (1973). The 
comparisons for 5% damping are shown in Figs. 4.31 and 4.32. The 
comparison is also summarized in Tables 4.5 to 4.8 for 2, 5 ,10, and 20 
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percent danping. 

Figure 4.31 for mean acceleraticsn and velocity anplification 
factors, clearly reveals that mean acceleration anplificaticn obtained 
from the Indian data is higher than that for data from earthquakes 
elsewhere. On the other hand, mean amplification in the velocity regicai 
is lower than those obtained by Nevmark-Hall and Mohraz. Also the 
velocity amplification factor reduces sharply with period in the 
Himalayan ground motions. However, mean-plus-one-standard-deviaticn 
acceleration amplification from Indian data is comparable (Fig. 4.32) 
with those obtained by Newmark-Hall, Mohraz and the one recoimended by 
NRC guide. This means that the standard deviation in acceleration 
amplification in the acceleratioi region is much lower for the Himalayan 
data than for data from elsewhere. 

4. 9 SPECTRUM INTENSm 

Peak ground acceleration is the most commonly used grovmd motion 
parameter in design. However, it has come under much criticism for its 
lack of correlation with observed structural performance. Hence, many 
other parameters have been proposed to measure the intensity of 
earthquake ground shaking. One such parameter is spectrum intensity, 
proposed by Housner (1959), Spectnan intensity (SI) is defined as the 
area under a pseudo-velocity spectnm curve between the periods 0.1 to 
2,5 sec, i.e. 


.5 

^ (T,^) dT 

Where S^ is PSV and P the percentage of critical damping. 



(4.6) 
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Tso et al. (1992) studied the correlation between A/V ratio and 
spectral intensity, with ^ being 5%. They observed that mean spectral 
intensity evaluated from PSV curves normalized to a comnon PGR of l.Og 
for low R/V (R/V less than 7.9 sec group of records is over four 
times higher than that for the his^ R/V grorq? (R/V greater than 11.8 
sec~^). Further, they show that SI correlates well with peak ground 
velocity (V), irrespective of R/V ratio (Fig. 4.33). They have obtained 
an empirical relationship between SI and V as 

SI = 3.3V (4.7) 

For the- data set used in this study, a similar correlation between 
SI and V is obtained. Spe<*rum intensity for 5% damping, for all the 
horizontal component records used in the stxjdy are presented in Rppendxx 
A for ready reference. Fig. 4.34 shows the plots of SI versus peak 
ground velocity, for the six earthquakes; also shown in the figure are 
the best-fit linear (in the form y = m x) relationships. Bnpirical 
relationships for eadi of the six events are presented below 
SI = 2.65V (For "/^ril”, M = 5.5) 

SI = 2.19V (For "Sept", M = 5.2) 

SI = 3.36V (For "May", M = 5.7) 

SI = 3.00V (For "Feb" M = 5.8) (4.8) 

SI = 3.56V (For "Aug", M =7.2) 

SI = 3.98V (For "Oct", M = 6.6) 

SI = 3.55V (For All records) 

The relationships indicate that, in general, slope of straight line 
correlating SI with V increases with increase in magnitude of 
earthquake. However, no definite relationship in slope and earthquake 
magnitude can be concluded from this study. 
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Table 4.1 Spectrum amplification factors for horizonta 
elastic response (Newmark and Hall 1969). 


Damping, 

% Critical 

One Sigma (84.1%) 

Median (50%) 

A 

[v _ 

D 

A 

V 

D 

0.5 

5.10 

KEOi 

3.04 

3.68 

2.59 

2.01 

1 


3.38 

2.73 

3.21 

2.31 

1.82 

2 

3.66 

2.92 

2.42 

2.74 

2.03 

1.63 

3 

3.24 

2.64 

2 . 24 

2.46 

1.86 

1.52 

5 

2.71 

2.30 

2.01 

2.12 

1.65 

1.39 

7 

2.36 

2.08 

1.85 

1.89 

*1.51 

1.29 

10 

1.99 


1.69 

1.64 


1.20 

20 

1.26 

1.37 

1.38 

1.17 

1.08 

1.01 


A : Acceleration 
V ; Velocity 
D : Displacement 


Table 4.2 Equations for Spectriim Amplification Factors 
for Horizontal Motion (Newmark and Hall 1969) 


Quantity 

Cumulative 
Probability, % 

Equation 

Acceleration 

Velocity 

Displacement 

84.1 

4.38- 1.04 InO) 

3.38- 0.67 InO) 
2.73-0.45 InO) 

Acceleration 

Velocity 

Displacement 

50.0 

3.21-0.68 InO) 
2.31-0.41 InO) 
1.82-0.27 InO) 


p : Fraction of critical damping 
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Table 4.3 Sutnnary of mean acceleration anplification factors 
for the six Himalayan earthquakes. 


Event 

period 

limits 

(sec) 

Percentage of critical danping 

1 0^5 1 

2.0 

5.0 

1 10.0 1 

l^ril 

0.13-0.24 

■B 

3.3 

2.6 

mgm 

Sept 

0.12-0.21 

B 

3.3 

2.6 

2.1 

May 

0.12-0.27 

5.8 

3.9 

2.8 

2.0 

Feb 

0.11-0.24 

mil 

3.7 

2.6 

mm 

Aug 

0.11-0.27 


3.6 

2.5 

1.9 

Oct 

0.14-0.29 

5.4 

3.8 

2. -7 

mm 

All 

records 

0.12-0.26 

mm 

3.5 

2.6 

2.0 


Table 4.4 Surmary of mean velocity anplification factors for 
the six Himalayan earthquakes. 
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Table 4.5 Catparitive sunrory of mean acceleration anplification 
factors. 


Percent of 

critical 

damping 

S^esent 

stxidy 

NeMnark & Hall 
(1969) 

MoBfkS (1976) 
for rock sites 

Period, 
0.12-0.26 sec 

Period, 
0.12-0.50 sec 

Period, 
0.12-0. 39 -sec 

2 

3.5 

2.74 

2.57 

5 

2.6 

2.12 

1.98 

10 

2.0 

1.64 

1.56 

20 

1.5 

1.17 

1.22 


Table 4.6 Conparitive sximary of mean velocity anplification 
factors. 


Percent of 

critical 

damping 

Present 

stiidy 

NeMonark & Hall 
(1969) 

Mcduraz (1976) 
for rock sites 

Period, 
0.26-2.60 sec 

Period, 
0.50-3.33 sec 

Period, 
0.33-3.33 sec 

2 

2. 6-1. 3 

2.03 

1.57 

5 

2.0-1. 1 

1.65 

1.28 

10 

1.4-0. 9 

1.37 

1.04 

20 

1.1-0. 9 

1.17 

0.81 
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Table 4.7 Ccnparitive sumary of mean-plus-caae-standard-deviation 
acceleration airplification factors. 


Percent of 

critical 

danpin? 

Present 

study 

Neunark & Hall 
(1969) 

Mdhraz (1976) 
for rock sites 

NRC 

Period, 
0.12-0.28 sec 

Period, 
0.12-0.50 sec 

Period, 
0.12-0.33 sec 

Period, 
0.11-0.40 sec 

2 

4.6 

3.66 

3.80 

3.54-4.25 

5 

2.9 

2.71 

2.82 

2.61-3.13 

10 

2.4 

1.99 

2.11 

2.27-2.72 

20 

1.7 

1.26 

i 

1.54 

1.90-2.28 


Table 4.8 Ccnparitive s\innary of mean-plxis-ane-standard-deviation 
velocity anplification factors. 


Percent of 

critical 

danping 

Present 

study 

NeHtrark & Hall 
(1969) 

Mohraz (1976) 
for rock sites 

Period, 
0.28-2.54 sec 

Period, 
0.50-3.33 sec 

Period, 
0.33-3.33 sec 

2 

3. 6-1. 9 

2.92 

2.44 

5 

2. 6-1. 6 

2.30 

1.90 

10 

1. 9-1.1 

1.84 

1.48 

20 


1.37 

' 1.11 





































Spectral Velocity, In/sec 



Fig. 4.1 SDOF System, subjected to earthquake ground motion 



Frequency, Hz 


Fig. 4.2 A typical response spectrum with a tripartite plot. 

Also, marked are the values of peak ground 
acceleration, peak ground velocity, and peak ground 
displacement (Gupta 1930). 






Idealized undamped velocity spectrum curves 
illustrating the effect of nagnitvide and distance, 
curve A, 25 miles from centre of large earthquake; 
curve B, 70 miles from centre of large shock; curve 
C, 8 miles from centre of small (M=5.3) shock 
(Housner 1970). 


Fig. 4.3 









Average Acceleration Coefficient (5a/g ) 
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PtRIOO. IN ifcCONDS 

0-5 0.2 0.1 0 05 0 02 0 01 



FREQUENCYJN HERTZ 


Fig. 4.6 "Standard" site independent horizontal response 

spectra modified from Newmark and Hall 1969 (Ilavn 
1980). 
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PERIOD, IN SECONDS 
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Fig. 4.7 


Horizontal design spectra and relative values of 
spectrum amplification factors for control points, 
adopted by NRC regulatory guide 1.6 (1973) (Hays 
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PERIOD. IN SECONDS 



FREQUENCY, IN HERTZ 


% of 

critical 

danping 

Accel erat ion 

Displacement 
D(0.25 Hz) 

A' (50 Hz) 

A(33 Hz) 

B(9 Hz) 

C(3.5 Hz) 

0.5 

0.67 

1.0 

4.96 

5.67 

2.17 

2.0 

0.67 

1.0 

3.54 

4.05 

1.67 

5.0 

0.67 

1.0 

2.61 

2.98 

1.37 

7.0 

0.67 

1.6 

2.27 

2.59 

1.25 

10.0 

0.67 

1.0 

1.90 

2.17 

1.13 


Vertical design spectra and relative values of 
spectnm airplification factors for control 
points, adopted by NRC regulatory guide 1.6 (1973) 
(Hays 1980). 


Fig . 4.8 





Fig. 4.9 


Average acceleration spectra for different site 
conditions (Seed et al. 1976). 
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Average 

acceleration amplification for 2 percent 

damping 

(Mohraz 1976). 
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Fig. 4.23 Mean-plus-one-standard-deviation spectra for the 
three earthquake magnitude (M) range. 



JMean"pl us“One“Standard "devi a t i on horizontal spectra 
considering all records together for 0, 0.5, 2, 5, 

10, and 20% damping. 


Fig. 4.24 



Period (sec) 

Pig. 4.25 OCV of acceleration spectra for 0, 5, aiKj 20% danping 
considering all horizontal records together. 



Coirparison of 18:1893-1984 spectra with mean spectra 
from six Himalayan earthquakes. 
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4.27 Cotparison of 5% danping spectral shapes obtained 
from data with that of 18:1893-1984 and UBC (1991). 
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Fig. 4.28 Ccrrparison of mean-plus-one-standard-deviation (m +a) 
spectra fran the Himalayan earthquakes data with 
"standard spectra" for rock sites reccrrrrended by Ml® 
c in onik damoina) . 




Fig. 4.30 


PSV spectra for 0, 2, 5, 10, zund 20% 

indicating straight line representation in 
region. 
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Fig. 4.31 


Conparison of mean airplification factors (a) 
acceleration anplification and (b) velocity 
anplification. 
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(a) 



ror™.;,rison of mean-plus-one-standard-deviation 

SifJcation factors (a) acceleration amplification 
and (b) velocity amplification. 


Fig . 4.32 
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Fig. 4.33 


Correlation between spectnm intensity and PGV for 
three groups of records having low, intermediate, and 
high A/V ratios. 
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Fig. 4.34 Correlation between spectrum intensity (SI) and PGV 
for the six Himalayan earthquakes. 












CHAPTER V 

DUCTILITY REDUCTION FACTOR 


5.1 INTROroCTIOK 

Elastic response spectrum is obtained by assuming that the 
structure behaves linearly elastic ev®a during severe earthguake ground 
motion. However, it will be impractical, or at least uneconomical, to 
design a structure to remain elastic during severe earthquake ground 
motion of rare occurence. H&ice, the seismic design specifications for 
structures must incorporate the inelastic response appropriately. It is 
well understood and recognized that when inelastic deformations are 
permitted in design by providing adequate ductility, the structure can 
be designed for a much lower seismic force than otherwise. 

The codes, particularly in the earlier years, provided for seismic 
design spectrum which was much lower than the elastic spectrum for the 
expected severe ground motion of rare occxirence. Such an approach 
implicitly assumed the inelastic response of structure. However, the 
current trend of seismic code develoEment is to prescribe a realistic 
level of elastic design spectrum correspcnding to severe ground motion. 
Force obtained by this spectrum is then reduced by a "response reduction 
factor" (R). 18:1893-1984 still has the first approach; many recent 
codes (ATC-1978, UBC-1991, NEHRP-1991, etc.) have adopted the later 
approach. The response reduction factor primarily consists of two 
factors. One factor is to reduce the elastic demand force to the level 
of maximum yield strength of the structure, this reduction basically 
depends on ductility of the structure and hence on its energy 
disRination capacity. This factor is called the ductility reduction 
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factor, R^. The other factor, called overstrength, is due to the 
overs tresngth inherently introduced m the cod€~designed structures 
is defined as the ratio between the inaxinun lateral strength of the 
structure and the code prescribed unfactored design base shear force. 

Fig. 5.1 shows idealized linearly elastic-perfectly plastic (EPP) 

global response of a building in tenns of base-shear coeffecient versus 

roof displacement relationship during severe ground shaking. Also shown 

in the figure is the response of building, if it remains elastic. The 

maximum base shear that develc^ in the structxare, if it were to remain 

elastic, is C ,W; where W is the total weight of building. During 

nonlinear response, the maximum roof displacement of the structure is 

A .If the global ductility demand H, defined as 
max 



can be achieved through appropriate design and detailing of structure, 

then such a building needs to be designed only for maxinun yield 

C 

strength of C^W, where Cy = . Hence, R^ plays an important role 

in evaluating the design force cn a structure. 

Ductility reduction factor (DRF, R^) depends on natural period and 
damping of the structure, its ductility, its load-displacement 
behaviour, and on the ground motion characteristics. In the present 
study, the structure has been assumed as a SDOF elastic-perfectly 
plastic (EPP) system. The value of DRF has been taken as the ratio of 
maxinun elastic response (PSA) to the maximum inelastic response (PSA) 
with ductility M for a SDOF EPP oscillator with 5% damping, i.e. , 
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%(T!rU) 


PSA (T, C = 5%)^^^ 
PSA (T, C = 5%, ^i) 


(5.2) 


In this chapter, the value of KCF has been studied for the 
Himalayan ground motion records, and compared with studies available in 
the literature for ground motions recorded elsewhere. 

5.2 REVIEW OF THE PAST STIJDIES 

The first attempt to relate R^^ with M was made by Newntark and Hall 
(1973, 1982) for SDOF systems with EPP force-deformation model. Based on 
analytical studies, they concluded that: 

- In structures of natural period less than 0.1 sec, ductility does not 
help in reducing the response. .Hence, for such structures, no 
ductility reduction factor should be used. 

- For moderate period structures, corresponding to acceleration regim 
of elastic response spectrum (T = 0.1 to 0.5 sec), the energy absorbed 
by an inelastic structure at its maximun displacement approximates 
that absorbed by an elastic system. This results in the ductility 
reduction factor as (Fig. 5. 2a) 

R ^= m - 1)^/2 ( 53 ) 

- For relatively long period structures, corresponding to velocity 
region (T = 0.5 to 3 sec) and displacemait region ( T = 3 to 10 sec) 
of elastic response spectrum, the maximum displacement of elastic and 
inelastic systems is nearly equal. This gives the value of ductility 
reduction factor as (Fig. 5. 2b) 

R^ = M 


(5.4) 



Figure 5.3 shows the DRF spectra proposed by Newmark and Hall (1973). 

Riddell and Newmark (1979, 1980) conducted a statistical analysis 
of the response of SDOF systems with elasto-plastic, bilinear, and 
stiffness degrading load-deformation model. For EPP, SDCF systems they 
obtained empirical formulae for R^ depending on ductility (m) and 
viscous danping ratio (C) for different period ranges. They observed 
that: 

- The ordinates of the mean inelastic spectra do not vary significantly 
whan various nonlinear models are used; differences occur mainly for 
intermediate frequencies and large ductility, and are practically 
negligible at the low and high frequency ends of the spectnin. 

- Use of the elasto-plastic idealization provides, in almost every case, 
a conservative estimate of the average response to a nvmber of 
earthquake motions. 

- It is particularly significant that, on the average, the stiffness 
degradation phenomenon is not as critical as one might expect a 
priori . 

Using ten "standard" ground motions from five differait seismic 
events, Mahin and Bertero (1981) have computed displacemKit ductility 
demands for EPP SDOF systems designed using Newmark-Hall inelastic 
design response spectra with 5% danping. They excluded that 
displacement ductility denand on the average is less than the values 
specified by Newmark-Hall method. However, Newmark-Hall method becomes 
less conservative as the values of the specified ductility and viscous 
damping are increased. 

Briseghella et al. (1982) have proposed a relation betweai and R^ 
as a function of period based on a study of SDOF systans, for Takeda 
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model of load-deformation relationship. Tliey have used artificial time 
histories conpatible with the power spectral density of a recorded 
motion during the 1976 Friuli earthquake. 

For EPP load-deformation relationship, Lai and Biggs (1980) have 
investigated the sources of variability of inelastic response spectra, 
due to strong motion duration, ductility, and danping. From time history 
analyses, they concluded that the inelastic response spectra are not 
significantly dependent on strong ground motion duraticm. However, they 
added that this conclusion is valid only when the strong ground motion 
with varying durations are coripatible with the same prescribed elastic 
response spectnm. They also found that the Newmark-Hall procedure for 
predicting the inelastic response of a 5% damped EPP system is 
unconservative, on the other hand for 2% damping the Newmark-Hall 
approach predicts a conservative inelastic respcmse. Fig. 5.3 shows the 
DRF spectra proposed by them. 

Pal et al . (1987) have developed yield displacement spectra, 
constant strength spectra, constant displacement spectra, reduction 
factor spectra, inelastic acceleration spectra, and inelastic yield 
displacanent spectra for three artificial and two actual earthquake 
time-histories using EPP, bilinear, and stiffness degrading models. They 
concluded that the design seismic loads for strong motion earthquake can 
be reduced below those required for elastic analysis by a factor of 
about 1.1 to 15 depending r^on the time period of the structure and the 
desired displacement ductility. Further, the reduction factors due to 
stiffness degrading model are only slightly higher than those for EPP 

nodels. This can be seen in Fig. 5.5. 

Miyann et al. (1988) have calculated the inelastic response of 
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SDOF systems for bilinear, and ’’origin" oriented models of 
load-deformation relationship. The study is based aa artificial 
earthquake ground motions of which the target spectmn is elastic 
response spectrum. They have expressed as functim of ductility 
factor, natural frequency of system, and shape of the elastic response 
spectra. 

Riddell et al. (1989) studied non-linear behaviour of EPP SDOF 
system, particularly for short period range, using 53 ground jnotion 
records, and concluded that reduction factor tends to a value very close 
to unity for very small periods and to a constant value larger than one 
for periods over 0.4 sec. The idealized bilinear relation between and 
T for different values of ductility is shown in Fig. 5.5. 

Using bilinear load-deformation relationship, Miranda (1992) has 
reported a statistical study on SDOF systems from 124 ground motion 
records from different types of site. The- study shows that local site 
conditions significantly affect the value of R^. Fig. 5.7 shows the 
variation of ductility reduction factor with time-period for rock sites 
and for soft soil sites. Notice that while the plot for rock sites is R^ 
versus T, for soft soil the plot is R^ versus T/T^. Here T^ is the 
predominant period of the site defined as the period correspcnding to 
the maximum spectral velocity. 

Vidic et al. (1992) state that peaks in the DRF-spectrum correspond 
to peaks in the elastic acceleration spectrum and that, for this reason, 
smooth DRF-spectra are applicable in combination with smooth elastic 
design spectra. Based cxi the response of an oscillator with Q-hysteretic 
model with 5% daitping, they show the influsice of gromd motions (Fig. 

5.8). 
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Thus, several studies (eg., Riddell et al. 1989, Vidic et al. 1992) 
canclvde that the value of R^ 

(a) in the short period range depends on T and m; it increases linearly 

with T, starting with M = 1 and going v 5 )to and 

(b) in the medium and long period range, depends on U but is independent 
of T; its value is around 


5.3 METHOD CP MIRLYSIS 

IHF spectra for all the horizontal component records were obtained 
for SDOF, EPP systems using Newmark-beta method. The displacem^t 
ductility value chosen was 2, 3, 5, and 10. The computations were 
carried out for 41 values of natural periods (0.001, 0.055, 0.01, 0.02, 

..., 0.09, 0 . 1 , 0 , 2 , .../ 3 * 0 )- 

mean-minus-one-standard-deviation, and coeffecient of variation (C30V) of 
ductility reduction factors for individual events and for the average of 
six events have been evaluated. 

5.4 RESULTS 

Mean CRP spectra in the six events for displacement ductility of 2, 
3, 5, and 10 are shown in Fig. 5.9. It is seen that, irrespective of the 
value of ductility, for very small periods, DRF is close to unity. 
However, for periods upto 0.2-0.3 sec, the spectral values build up and 
later remain fairly constant upto period of about 2.7 sec. R slight 
increase in PRF value is observed in the period range 2.7 to 3.0 sec. 


5.4.1 Mean DRF Spectra 





The mean DE^F spectra for different evaits are cotipared in Fig. 5.10 
for the four values of ductility. It is seen that for low ductility 
(less than 5) spectral shapes are quite consist^t for different events 
with only nominal variations. However, for higher value of ductility the 
variatiOTi in mean ERF spectrum tends to increase. This is in line with 
the observation by Riddell et al. (1989), that "the inelastic response 
of SDOF systems with large inelastic incursions are more affected by the 
characteristics of accelerograms than those with small inelastic 
incursims". 

The mean ERF spectra in "Sept" achieve an unusual peak at natural 
period of around 0.2 sec to 0.3 sec. However, beyond 0.4 sec the mean 
ERF spectra of this event are lower than thosefor all the other events. 
Notice that while the mean elastic spectra of this earthquake (Fig. 
4.16) shows a sharp drop at natural period of around 0.4 sec, the 
elastic response spectra of this earthquake do not show large 
amplifications in the period range of 0.2 sec to 0.4 sec. Hence, the DRF 
spectra of this evoit do not reflect the elastic response spectrum 
characteristics; this is in contrast with the observations of Vidic et 
al. (1992). The DRF spectra are the highest for event "April" as 
conpared to other events. Interestingly, the amplification factors in 
the acceleration and velocity regions of elastic spectra for this evmt 
are the lowest. 


5.4.2 OOV and Mean-^^inxis-One-Standard-Deviation Spectra 

The coeffecient of variation in ordinates of DRF spectra, for 
different ductility, for the six events is shown in Fig. 5.11. The 
spectra is shown for period range of 0.1 to 3.0 sec. It is vary clear 
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that the COV is higher for large value of ductility and vise versa. 
Also, irrespective of event, CX)V remains mostly between 0.2 to 0.4, for 
ductility value between 10 to 2. These observations can be clearly seen 
in Fig. 5.12, which shows the average of OOV spectra from each event. 
This has impotant implications for design of important structures where 
the ORF to be used is more like mean-minus-one-standard-deviation than 
the mean value. Such a consideration is not being made in the profession 
at present. 

Fig. 5.13 shows the of mean-ndnus-CHie-standard-deviation spectra in 
the six events. In general, shape of mean-ndnus-one-standard-deviation 
spectra is same as that of mean D3RF spectra. 


5.4.3 Idealisation of Mean and Mean-Minus-One-Standard-Deviation ^?ectra 
Fig. 5.14 and Fig. 5.15 show the mean and 
mean-minus-one-standard-deviation spectra obtained by averaging the 
corresponding spectra obtained in the six events. These shapes are 
idealised as bilinear curves and are shown in the same figure. Similar 
shapes were obtained by Riddell et al. (1989) (Fig. 5.6). The relation 
between DRF and period for idealised shapes of mean and 
mean-minus-one-standard-deviation spectra can be expressed as; 

* 


V ^ 

= 1.0 + ; — ^ 


for 0 i T < 0.2 sec 
for 0.2 ^ T < 3.0 sec 


(5.6) 


The values of R* for ductility 2, 3, 5, and 10 are given in Table 5.1. 
It can be seen from Fig. 5.9, that for longer periods (T > 2.7 sec), for 
high ductility (fi i 5), a slightly higher ductility reduction factorthan 



tut indicat^J by idealised shape is affordable. Thus, idealized sha^ 
in this period range is conservative. 

5.5 COMPARISON OP RESULTS 

Fig. 5.16 shows the idealized DRP spectra obtained from the data in 
conparisQii with the DRF spectra proposed by, Newmark and Hall (1973), 
Lai oind Biggs (1980) and Riddell et al. (1989) studies. The figure 
cleat ly inviicates that, the variation in different methods exists mostly 
in th*;' petiai tanye of 0.2 to about 0.6 sec. The spectra proposed by 
Lai -Biggs is ivitJ£;t con.servative. The Newmark-Hall spectra matches quite 
well with the pioposed spectra for, periods longer than 0.6 sec. But, 
for -ihiviter periods the Newmark-Hall spectra is quite lower than the 
dpiluined spectra. The proposed idealised shape of DRf spectra is similar 
to the shape proposed by Riddell et al. (1989). However, there are some 
differences in the period from which R^ becomes constant. For hig^r 
ductility (5 s m ^ 10), Riddell et al. suggests this period as 0.4 sec 
as agairist 0.2 sec proposed from the data. Due to this difference the 
slope of initial linear poiticai of the two bilinear curves vary. 



Table 5.1 


for mean and 

standard-deviation KIF spectra. 


Ductility 

For mean 
spectra 

For m-<y 
spectra 

2 

2.2 

1.7 

3 

3.2 

2.2 

5 

4.8 

3.2 

10 

8.7 

5.2 


m-c ; mean-minus-one-standard-deviation 
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Fig. 5.4 Ductility reduction factor spectra (Lai and Biggs 
1980). 




Coni>arison of ductility reduction factor spectra for 
elasto-plztstic and stiffness degrading models (Pal et 
al. 1987). 


Fig. 5.5 
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Pig. 5.7 Mean of strength reductions due to ncai-1 inear 
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Fig. 5.10 Cannparison of for each ductility for the six 
Hinalayan earthquakes. 
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CHAPTER VI 

SUMMARY AND CONCLUSIONS 


When formulating design criteria for aseismic design of structures, 
it is necessary to specify the ground motion that the structures must be 
designed to resist. The present knowledge in earthquake engineering and 
related sciences is not enough to know the exact characteristics of the 
future ground motion. However, data available from the past earthquakes 
at the site could be effectively used to estimate the earthquake ground 
motion in future earthquake. 

In the present study, an atteirpt has been made to study the data 
available from the six Himalayan earthquakes, which occurred during 
1986-1991. A data set consisting of 198 horizontal component and 99 
vertical component records is used to study the ground motion 
characteristics such as peak parameters, duration of strong ground 
motion, shape of response spectrum, and ductility reduction factor. The 
results obtained were compared with those available from similar studies 
on earthquakes elsewhere, and with currently adopted design practices. 

Following are the salient conclusions regarding characteristics of 
ground motions in the Himalayas from the six earthquakes recorded in the 

region: 

- The A/V ratio for the ground ' motion in the Himalayan region is 
unusually high. Hie mean Vv is nore like 20.0 sec’^ for rook sites; 
this is quite high when compared to the value of 12-15 sec reported 
for earthquakes elsewhere. The mean-minus-one-standard-deviation A/V 
value is obtained as 12 sec"^. 

- The variation in A/V ratio with the distance from shock-source is 



somewhat insignificant for the Himalayan earthquake ground motions. 
Thus, high A/V value is obtained even at large distances. This 
suggests that attenuation characteristics of peak ground motions in 
this region is different from those reported in literature. 

- The mean and mean-plus-one-standard-deviation value of AD/V^ are 

obtained as 4.6 and 8.0, respectively. No clear trend in variation of 
2 

AD/V with distance is observed. 

- The mean and mean-plus-one-standard-deviation A^Ajj value is obtained 
as 0.51 and 0.70, respectively. A conservative value of 0.67 or 0.75 
is suggested as against the currently adopted value of 0.5 in Indian 
code. 

- A/V is found to be good indicator of frequency content of ground 
motion and hence the relative shape of spectra. 

- Duration of strong shaking, as evaluated by Trifunac and Brady (1975) 
definition, shows no definite relationship with A/V or with distance 
from shock-sotirce, though sxich relationships are reported in 
literature. 

- Significant variation in elastic spectral shapes in intermediate 
period range are observed depending on the iragnitude of earthquake. 
Mean elastic spectral shapes for three ranges of earthquake 
magnitude; (1) small magnitude (M < 5.5), (2) moderate magnitude (5.5 
s£ M ^ 6.5), and (3) large magnitude (M > 6.5) are obtained. This could 
be helpful in deciding the shape of design spectra. 

- Spectral shapes do not vary much with distance for moderate magnitvde 

A. 

earthquakes within moderate epicentral distances. 

- Both, the ratios of peak parameters and the elastic response spectra 
indicate that with distance the attenuation of waves of different 
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fre<juencies is not so different as elsewhere. Hence, for projects in 
Himalayas, the enpirical relationship for peak parameters with 
magnitude, epicentral distance, etc., of elsewhere are not applicable. 
The design spectral shape provided in OBC (1991) is quite conservative 
if applied to Himalayan earthquake ground motion as it completely 
envelops the obtained mean shape. The shape of design spectra 
currently adopted in 13:1893-1984 is quite different from the mean 
shape obtained from data. The code spectra underestimates the spectral 
ordinates in period range of 0.1 to 0.35 sec and overestimates in the 
period range greater than 0.35 sec. Revision of currently adopted 
shape in Indian code is urgently required. In order to bring the mean 
spectral shape of Indian code at par with the shape obtained from the 
data, the acceleration amplification in the period range of 0.1 to 
0.35 sec may be increased, while the presently used shape for longer 
periods is conservative and can be kept as it is. 

Shapes of AERB "standard spectra" for rock sites and that of 
mean-plus~one-standard-deviation spectra from the data are in good 
agreement. This indicates that spectral shapes suggested by AERB are 
quite satisfactory for the structures which are usually designed using 
mean-plus-one-standard-deviation spectra . 

The acceleration amplification region mostly ranges between natural 
periods of 0.12-0. 26 sec, while velocity region ranges between periods 
0,26 to 2.6 sec. Mean acceleration amplification factor for the data 
is higher than that for data from earthquakes elsewhere. On the other 
hand, mean airplification factor in velocity region is lower. However, 
mean-plus-one-standard-deviation acceleration amplification from 
Indian data is comparable with those obtained in other studies. 
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Straight line correlation between spectmn intensity and peak ground 
velocity are established. It is found that the slope of straight line 
correlating spectrun intensity with peak ground velocity increases 
with the increase in earthquake magixitude. 

A relationship between ductility reduction factor and time period of 
elasto-plastic single degree of freedom structure is proposed. Except 
for some difference for periods less than 0.4 sec the proposed spectra 
is quite similar to the one proposed by Riddell et al. (1989). 

Based on the preset work, following area for research in future are 

identified. 

- A/V is an excellent indicator of frequency content. Hence, design 
recanmendation can be based on VV ratio. This approach is adopted by 
National Building Code of Canada (NBCX; 1990). Hence, study of A/V 
ratio to work out design reccmnendations can be carried out. 

- Amplification factors in displacement region could not be evaluated 
because of the unexpected shape of spectra in this region. A detailed 
investigation on this is needed. 

- The present work on ductility reduction factor is limited to 
elasto-plastic single degree of freedom systems, this work could be 
extended for other load-deformation models and for multi degree of 
freedom systems. 

- Now better strong motion instrumentation is available, to make full 
use of it the ground motion parameters like magnitude, location of 
epicentre, etc., are needed with better precision. This will be very 
helpful in developing the attenuation characteristics of the peak 
ground motion parameters. 

- In India, detailed studies on analysis of available strong motion data 



are only few, so many other parameters and characteristics could be 
studied so that one could reach to the current state-of-art in 
Califorma or Japan. 
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APPENDIX A 


Table A-1 For April 26, 1986 earthquake recorded in Kanqra Array. 


SR. 

NO. 

LOCATION 

'“epicentral 

DISTANCE 

(Km) 

dirn! 

PEAK 

GROUND 

ACCL, 

(q) 

PEAK 

GROUND 

VEL. 

( mm/s ) 

PEAK 

GROUND 

DI5PL. 

(mm) 

DURATION 
TRIFUNAC 
?x BRADY 
(1975) 

(sec ) 

SPECTRUN 
INTENSITY 
FOR 5% 
DAMPING 

(mm) 

1 

BHANDLKHAS 

A1 


14.8 

L 

0.145 

83.13 

20.08 

2.82 

174.92 



A2 


22.3 

T 

0.125 

90.75 

5.85 

3.16 

147.23 



A3 

= 

26.3 

V 

0.020 

19.03 

6.33 

9.04 

61.93 



A4 

» 

24.9 


‘ 





2 

BAROH 

A1 

n 

10.1 

L 

0.058 

36.00 

6.55 

4,34 

G9.03 



A2 

a 

13.1 

T 

0.057 

27.00 

4.58 

6.74 

101.74 



A3 

a 

22.2 

V 

0,020 

15.07 

3.64 

10.50 

57.85 



A4 

= 

23.2 







3 

BHAWARNA 

A1 

= 

14.7 

L 

0.036 

11.70 

2.73 

9.04 

58.39 



A2 

a 

19.1 

T 

0.035 

19.21 

4.20 

8,08 

75.72 



A3 

= 

26.3 

V 

0.036 

9.34 

21.13 

5.98 

36,85 



A4 

a 

25.5 







4 

DHARANSHALA 

A1 

a 

9.9 

L 

0.176 

72.97 

7.76 

2.60 

188.37 



A2 

a 

11,3 

T 

0.186 

94.90 

24.79 

3.94 

379.46 



A3 

a 

5.0 

V 

0.0B3 

27.39 

4.19 

4.34 

99.61 



A4 

a 

2.8 







5 

JAWALI 

A1 

a 

36.5 

L 

0.015 

19.20 

6.90 

: 15.92 

78.62 



A2 

a 

29.0 

T 

0,017 

15.18 

9.56 

13.99 

58.34 



A3 

r 

26.4 

V 

0.011 

14.39 

7.53 

16.14 i 

71.01 



A4 

a 

28.6 







6 

KANGRA 

A1 

a 

13.5 

L 

0.14B 

50,57 1 

5.23 

4.84 

143.84 



A2 

3 

5.2 

T 

0.111 

95.7''i 


3.?4 




A3 

a 

11.0 1 

V 

0.072 

31.93 

5.37 

5.40 

104.75 



A4 

a 

12.5 







7 

NABROTABAGWAN 

A1 

8 

3 . 6 

L 

0.149 

94.17 

13.09 

4.20 

220.38 



A2 

a 

6.4 

T 

0.080 

25.40 

6,14 ' 

5,14 

78.55 



A3 

a 

12.9 

1 

V 

0.050 

18.19 

6,51 

14.82 

89.39 



A4 

a 

12.5 







8 

SHAHPUR 

A1 

= 

20.4 

L 


59.21 

7.26 

2.70 

193.59 



A2 

a 

15.8 

T 

0.250 

147.80 

10.85 

2,00 

309.97 



A3 

a 

6.8 

V 

0.065 

28.04 

5-25 

6.40 

110.13 



A4 

a 

10.3 







9 

SIHUNTA 

A1 

a 

34.3 

L 

0.050 

26.72 

4.71 

6.20 

96.32 



A2 

a 

30.6 

T 

0.036 

33.88 

3,67 

6 . 66 

105.27 



A3 

a 

22.4 

V 

0.039 

27.73 

4.74 

9.34 

78.79 



A4 

3 

23.4 








# L • Longitudinal, T : Transverse, 


V : Vertical 



134 


T.bl. A-2 For Bfottmbtr 10, 1906 ..rthqu.k. recorded in Shillong arrty. 


SR- 

NOe 

LOCATION 

DISTANCE 

(Km) 

dirnT 

PEAK 

BROUND 

ACCL. 

<9) 

PEAK 

GROUND 

VEL. 

(mm/s) 

PEAK 

GROUND 

DISPL. 

(mm ) 

DURAfiON 

TRIFUNAC 

BRADY 

(1975) 

(sec ) 

spectrUn 

INTENSITY 
FOR 5% 
DAMPING 

(mm) 

1 

BAITHLANBSO 

A1 

u 

139.4 

L 

0.045 

20.58 

3.55 

10.38 

60.31 



A2 

s 

102.1 

T 

0.042 

12.07 

2.73 

10.78 

55.40 



A3 

= 

101.2 

V 

0.025 

10.91 

4.74 

11.86 

39.24 

2 

DAUKI 

A1 

a 

42.7 

L 

0.089 

32.40 

3.00 

3.92 

64.50 



A2 

a 

32.5 

T 

0.090 

37.33 

2.56 

5.30 

67.49 



A3 

s 

48.1 

V 

0,032 

14.14 

5.92 

14.26 

60.31 

3 

KHLIEHRIAT 

A1 

a 

ei.6 

L 

0.031 

11.30 

2.00 

10.56 

35.00 



A2 

= 

56.9 

T 

0.046 

18.41 

3.86 

9.60 

49.81 



A3 

s: 

54.4 

V 

0.001 

7.71 

1.95 

11.08 

37.01 

4 

NONBKHAW 

A1 

c 

65.4 

L 

0.055 

33.09 

9.32 

17.24 

103.41 



A2 

= 

35.8 

T 

0.093 

48.08 

5.47 

11.76 

120.43 



A3 

= 

28.3 

V 

0.034 

14.11 

6.83 

23.50 

64.77 

5 

NONGPOH 

A1 

a 

93.5 

L 

0.054 

21.13 

4.26 

8.32 

39.04 



A2 

s: 

56.4 

T 

0.056 

11.32 

4.59 

12,86 

35.94 



A3 

r 

33.8 

V 

0.033 

11.25 

2.04 

11.19 

27.90 

6 

NONGSTOIN 

A1 

s 

55.8 

L 

0.019 

10.04 

4.90 

7.76 

32.69 



A2 

w 

54.5 

T 

0.014 

6.40 

1.66 

7.54 

74.30 



A3 

a 

64.3 

V 

0.008 

5.95 

1.99 

7.82 

22.68 

7 • 

RANIHUR 

A1 

a 

134.9 

L 

0.040 

9.37 

2.02 

8.94 

27.08 



A2 

= 

103.8 

T 

0.048 

21.15 

2.65 

5,78 

39.76 



A3 

a 

90.2 

V 

0.023 

5.61 

1.40 

9.42 

27.87 

8 

PYNURSLA 

A1 

s 

37,4 

L 

0.093 

26.75 

3.24 

5.14 

83.82 



A2 

a 

15.0 

T 

0.076 

20.37 

4.55 

8.68 

62.32 



A3 

= 

33.6 

V 

0.030 

9.28 

3.80 

14.78 

35.17 

9 

SAI TSANA 

A1 

X 

104.4 

L 

0.110 

36.69 

4.90 

6.12 

74.97 



A2 

X 

68.3 

T 

0.140 

50.73 

4.92 

2.44 

106.37 



A3 

= 

50.1 

1 

0.062 

20.06 

3.22 1 

7.32 

54.22 

10 

UNNULONG 

A1 

X 

70.9 

L 

0.110 

26.53 

2.38 1 

5.82 

46.13 



A2 

a 

36.7 

T 

0.063 

12.10 

1.56 

9.56 

34.91 



A3 

m 

27.4 

V 1 

0.050 

14.25 

3.21 

11.42 

47.99 

11 

UNRONGSO 

A1 

X 

111.1 

L 

0,027 

9.68 

1.30 

8.26 

32.03 



A2 

X 

82.4 

T 

0.032 

11.84 

4.41 

9.82 

35.96 



A3 

X 

72.5 

V 

0.014 1 

1 

7.81 

2.26 

10.66 

32.72 

12 

UNSNING i 

A1 

X 

76.4 

L 

0.100 

28.72 

3.42 

4.26 

57.56 



A2 

X 

38.3 

T 

0.076 

25.67 

3.93 

8.78 

64 . 84 



A3 

X 

15.0 

V 

0.050 

10.88 

3.75 

14.50 

41.64 


# L « Longitudinel, T i TransverBe. V t Vertical 
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Fable A--'3 For May 


SR. 

NO. 


LOCATION 


BAITHLANGSO 


BANUN6A0 


BFRLONGFER 


1 \ 


4 


5 


A 


7 


8 


9 


10 


1 1 


12 


BOKAJAN 

DIPHU 

QUHmm 

HAFLONG 

HAJADISA 

HATIKHALI 

jtAISONG 

NONGPOH 

PANINUR 


18, 1987 earthquake recorded in Shillong array. 


EPICENTRAL 

DISTANCE 

(Km) 

DIRwt 

PEAK 

GROUND 

ACCL. 

<9> 

PEAK 

GROUND 

VEL. 

(mm/ s ) 

PEAK 
GROUND 
DISPL . 

(mm ) 

DURATION 

TRIFUNAC 

BRADY 

(1975) 

(sec ) 

SPECTRUM 
INTENSITY 
FOR 5*/. 
DAMPING 

< mm ) 

A1 

ss 

121.8 

L 

0.034 

15.53 

4,15 

16.84 

55.89 

A2 

= 

177.5 

T 

0.027 

21.56 

8.11 

18,50 

70.02 

A3 

- 

112.8 

V 

0.020 

12.85 

3.80 

18.14 

59.79 

A1 

55 

82,1 

L 

0,020 

12.22 

4.24 

22.58 

44.31 

A2 

3 

137.7 

T 

0.020 

22.00 

6.78 

23.48 

57.87 

A3 


73.8 

V 

0.019 

13.62 

3.32 

25.18 

65.10 

A1 

3 

54.0 

L 

0.072 

35,95 

5.51 

26.52 

110.19 

A2 

3 

109.5 

T 

0.090 

36 . 20 

8.95 

22.78 1 

128.03 > 

1*^ 

3 

46 , i:> 

V 

Q.mh 

19.24 

7.47 

30.62 

72.52 

A1 

3 

56,4 

L 

0.030 

25.96 

9.34 

18.88 

75.47 

A2 

3 

91.5 

T 

0.066 

34.28 

10.31 

12.04 

107.65 

A3 

= 

58.4 

V 

0.020 

13.20 

5.25 

21.92 

44,61 

A1 

ss 

53.7 

L 

0.086 

27.69 

5.80 

27.06 

94.39 

A2 

3 

104,7 

T 

0.073 

23.92 

11.54 

27.18 

90.79 

A3 

= 

49.6 

V 

0.054 

15.31 

6.65 

26.68 

78.30 

At 

3 

73.5 

L 

0.040 

23.04 

5.55 

12.46 

84.34 

A2 

3 

120.0 

T 

0,049! 

28.56 

5.51 

11.60 

90,31 

A3 

= 

64.0 

V 

0,019 i 

10.83 

3.38 

12.08 

53.12 

A1 

3 

77,9 

L 

0.055 

36.93 

i 

8.35 i 

6 • 62 

120.59 

A2 


119.0 

T 

0.035 

20.61 

5.38 

9,82 

107.89 

A3 

3 

69,3 i 

V 

0.015 

12.15 

3.06 

11.06 

62.22 

A1 


42.7 

L 

0.078 

38,07 

7.50 

12,16 

111,86 

A2 

- 

90.5 

T 

0.086 

28.07 

5.00 

11.24 

99.50 

A3 

3 

33.6 

V 

0.028 

15.34 

3.49 

13.40 

74.82 

A1 


62.3 

L 

0.031 

14.55 

5,55 

27.06 

60 . 68 

A2 


117.5 

T 

0.038 

17.71 

3.48 

27.26 

70.48 

A3 

= 

52.7 

V 

0.027 

18-58 

3,64 

28,86 

74.64 

A1 


52.2 

L 

0.042 

32.25 

7.84 

11,28 

103.30 

A2 

3 

90.4 

T 

0.061 

23.69 

9.42 

9.68 

109. 10 

A3 

3 

45.0 

V 

0.019 

13,54 

6.20 

13,68 

67.30 

A1 


188.7 

L 

0.017 

13.12 

4.88 

' 16.24 

35.71 

A2 

3 

243.4 

T 

0,047 

12.58 

5.93 

16.46 

51.92 

A3 

= 

179.0 

V 

0.013 

9.55 

3.13 

16.62 

54.64 

A1 


91 ,6 

L 

0.040 

18,00 

4, 10 

9. 12 

53.32 

A2 

3 

146.0 

T 

0.047 

18.11 

2.99 

9.92 

47.74 

A3 

= 

81.7 

V 

0.017 

16.79 

5.78 

9.52 

37,23 
. con td . 
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5R. 

siO. 

LOCATION 

EP I CENTRAL 
DISTANCE 

(Km) 

dirn!" 

PEAK 

GROUND 

ACCL. 

<9> 

PEAK 

GROUND 

VEL. 

(mm/s > 

PEAK 

GROUND 

DISPL. 

( mm ) 

’immrnr 

TRIFUNAC 

BRADY 

(1975) 

(sec ) 

SPECTRUfr” 
INTENSITY 
FOR 5% 

DAMPING 

(mm ) 

13 

SAITSAMA 

A1 

= 

133.9 

L 

0.037 

13.92 

3.36 

20.36 

39. 16 



A2 

= 

188. 1 

T 

0.050 

22.84 

8.91 

20.38 

80, 19 



A3 

s 

124.1 

V 

0.019 

9.29 i 

4.00 

21 .14 

45,12 

14 

UMRONGSO 

A1 

S 

106.3 

L 

0.021 

13.79 

2.76 

10.94 i 

44,37 



A2 

=: 

160.2 

T 

0.026 1 

14.50 

5.69 

9.90 

51.81 



A3 


98.3 i 

V 

0.016 

9.13 

2.91 

10.04 

41.11 


L : Longitudinal, T : Transverse, V ; Vertical 
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lb 1 < 

? A*-4 For Febn 

uary 

6 

, 1987 

earthquake recoded in Shi Hone Arrav. 



tPlCENTRAL 





DURATION 

SPECTRUN 

>R . 

L OC AT I ON 

DISTANCE 

dirn* 

PEAK 

PEAK 

PEAK : 

TRIFUNAC 

INTENSITY 

iO « 






GROUND 

GROUND 

GROUND : 

& BRADY 

FOR 57. 







ACCL. 

VEL. 

DISPL. : 

<1975) 

DAMPING 




(Km) 


(g) 

(mm/s ) , 

<mm) 

<sec ) 

( mm > 


BA I GAO 

A1 

= 

190,3 

L 

0.022 

14.14 

3.90 

11.28 

1 

56.19 



A2 

s; 

159,0 

T 

0.024 

20.64 1 

6.90 

10.72 

115.95 



A3 

s 

150.4 

V 

0.009 

11.271 

3.26 

10-58 

73.46 

> 

BAITHLANGSO 

A1 

s 

193.2 

L 

0.031 

13.65 

3.81 

10.06 

59,99 



A2 

= 

184.0 

T 

0.022 

12.26 

3.33 

9.94 

70.49 



A3 

- 

135.5 

V 

0.014: 

7.44 

1.76 

10.46 

29 . 40 


BANUNGAO 

A1 


223,1 

L 

0.016! 

5.08 

1 .30 

7.24 

27,01 



A2 


203.9 

T 

0 . 013 ! 

10.35 

2.30 

7.72 

; 33 . 88 



A3 


175.5 

V 

O.OIO! 

6.23 

2.33 

7.20 

26.98 

i- 

DAUKI 

A1 


103.7 

L 

0.027 

18.01 

2.50 

7.50 

55.74 



A2 


79,0 

T 

0.038 

18,37 

4.41 

6.78 

57.38 



A3 

= 

75.5 

V 

0.010 

8,57 

3.41 

8.62 

33 . 73 

) i 

GUNJUNG 

A1 

=s 

202.8 

L i 

0.036 

29,05 

6.26 

8.68 

1 103.77 

1 

j 


A2 


166.7 

T 

0-037 

24.44 

3.87 

7.96 

74.97 



A3 

= 

166,5 

V 

0.020 

12.87 

3.30 

10.78 

44.54 

> 

HAFLONG 1 

A1 

s: 

202.5 

L 

0.034 

19.42 

3.13 

8.16 

1 67,35 



A2 

a 

162.3 

T 

0.027 

14.02 

3.34 

8.04 

j 71,51 



A3 

a 

170.8 

0 

0.008 

7.22 

1 .99 

1 

1 

9.04 

1 22 . 00 

1 

r 

HATIKHALI 

A1 

a 

222,3 

L 

0.024 

8.78 

3.18 

8.98 

33.11 



A2 

a 

195.3 

1 T 

0.025 

9.87 

3.84 

10.18 

33.50 



A3 

= 

175.8 

1 

0.021 

7,29 

1 .98 

10,34 

28.00 


|KATAKHAL 

: A1 

a 

164,6 

L 

1 0.009 

1 10,64 

3,04 

8.82 

64.70 



1 A2 

a 

114.9 

T 

1 0.008 

j 7,75 

1 .60 

9.28 

32.11 



A3 


149,4 

V 

: 0.010 

j 

9,88 

4.00 

9.30 

44.81 

} 

KHLIEHRIAT 

A1 

a 

142,2 

L 

1 0.080 

47,46 

8.45 

3.90 

133.52 



A2 

a 

116,9 

T 

0.066 

27.87 

3,44 

6.12 

85.16 



A3 

= 

103.0 

V 

0,029 

1 15.91 

2,94 

9.30 

49.40 

,0 

MAWPHLANG 

A1 

a 

91 .7 

L 

0,081 

, 43,85 

7.97 

19.90 

98,18 



A2 

a 

93.5 

T 

0.059 

24.89 

10,14 

20-28 

100.55 



A3 

« 

42.0 

V 

0,036 

12,33 

4,86 

1 

23.74 

42. 19 

. 1 

NONGKHLAW 

A1 

a 

99.3 

L 

0.110 

42.94 

9.58 

24.20 

89,75 



A2 

a 

116.5 

T 

0,114 

52.63 

11.41 

13,52 1 

156.09 



A3 

= 

35.0 

V 

0,104 

1 35.10 

11.23 

16.58 

84,15 

l2 

NONGPOH 

A1 

a 

132.9 

L 

0.027 

15.89 

4,64 

12.28 

46.50 



A2 

s: 

144.7 

T 

0.086 

1 26.45 

1 3.75 

8.40 

74.80 



A3 

a 

68.4 

V 

0.038 

10.64 

1 4.02 

12.68 

40, 16 


. . con td . 
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SR. 

NO. 

LOCATION 

EP I CENTRAL 
DISTANCE 

(Km) 

DIRNt 

PEAK 

BROUND 

ACCL. 

(9> 

PEAK 

BROUND 

VEL. 

(mm/s ) 

PEAK 

BROUND 

DISPL- 

(mm ) 

DURATION 

TRIFUNAC 

Zi BRADY 
(1975) 

(sec ) 

SPECTRUM 
INTENSITY 
FOR 5*/. 
DAMPING 

(mm ) 

IB 

PYNURSLA 

A1 

r: 

96.1 

■B 

0.050 

21 .24 

4.99 

21.96 

72.09 



A2 

5S 

83.1 


0.031 

11.95 

4.93 

22.54 

51.90 



A3 

= 

59.4 

H 

0.015 

9.43 

3.63 

28.44 

47.79 

14 

SAITSAHA 

A1 

£ 

159.9 


0.066 

22.38 

4.50 

9.28 

71 .83 



A2 

£ 

148.6 

T 

0.058 

32.43 

4.37 

9.66 

75.86 



A3 

£ 

106. 1 

V 

0.032 

9.90 

2.37 

12.52 

48.50 

15 

SH It LONS 

A1 

£ 

109.1 

L 

0.047 

16.06 

5.16 

8.92 

57.06 



A2 

£ 

108-4 

T 

0.036 

15.39 

7.86 

9.00 

64.78 



A3 

= 

55.3 

V 

0.013 

14.42 

5,19 

10.64 

57.49 

16 

1 

UNhULONG I 

A1 

£ 

127.9 

L 

0.056: 

22.481 

5.69 

19.82 ! 

66 • 66 



A2 

£ 

115.3 

T 

0.055 

21.42 

6.35 

19.18 

81.74 



A3 

= 

80.0 

V 

0.024 

10.89 

4 . 37 ; 

22.30 

40.30 

17 

UHRONGSO 

A1 

£ 

171 . 1 

L 

0.046 

31.43 

3.06 

7.80 

81 .06 



A2 

£ 

146.7 

T 

0.037 

26.58 

6.67 

9.32 

95.65 



A3 


126.9 

0 

0.022 

12.26 

2.38 

11.20 

47.64 

18 

UHSNING 

A1 

£ 

205.5 

L 

0.040 

16-62 

3.14 

18.19 

51,45 



A2 

£ 

184.0 

T 

0.061 

36.93 

5.03 

14.54 

78,16 



A3 

£ 

155.5 

V 

0.018 

11.95 

3.55 

20.84 

53.25 


L : Longitudinal, T ; Transverse, V : Vertical 
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able A 5 For August 6, 1988 earthquake recorded in Shillong array 


SR. 

NO. 

LOCATION 

EP I CENTRAL 
DISTANCE 

(Km) 

dirnT 

PEAK 

GROUND 

ACCL. 

(g) 

PEAK 

GROUND 

VEL. 

(mm/s) 

PEAK 

GROUND 

DISPL. 

(mm ) 

DURATION 
TRIFUNAC 
&c BRADY 
(1975) 

(sec ) 

SPECTRUM 
INTENSITY 
FOR 5% 
DAMPING 

(mm ) 

1 

BAIGAO 

A1 

Si 

227.5 

L 

0.221 

65.68 

10. 19 

31 .919 

169.00 



A2 

s: 

226.7 

T 

0.144 

63.70 

9.80 

33.899 

182,74 



A3 

HZ 

165.0 

V 

0.053 

20.11 

5.85 

38. 180 

83.00 

2 

BAITHLAN6S0 

A1 

= 

268.2 

L 

0.154 

75.52 

12.36 

28,079 

256.94 



A2 

s 

266.2 

T 

0.165 

117,05 

15.45 

28,839 

420 . 49 



A3 

- 

200.3 

V 

0.082 

36. 12 

6.98 

46.099 

136.21 

3 

BANUN6A0 

A1 


227.6 

L 

0.093 

64.03 

8.96 

24,319 

203.34 



A2 


225.5 

T 

0.069 

49.80 

10.23 

22.979 

200,59 



A3 

= 

159.5 

V 

0.065 

25.06 

4.23 

23.359 

85.14 

4 

BERLONGFOR 

A1 


199.5 

L 

0.300 

217.24 

33.40 

21.960 

898 . 1 1 



A2 

= 

197.5 

T 

0.344 

228.19 

36.30 

20 . 000 

806.13 



A3 

- 

131.6 

V 

0.174 

90.26 

13.18 

41.599 

247.96 

5 

BOKAJAN 

A1 

3= 

167.2 

L 

0.151 

86.75 

19.91 

20,699 

345.83 



A2 

S 

163.8 

T 

0.224 

121.36 

20.41 

20 . 499 

369.27 



A3 

33 

99.5 

V 

0.148 

32.49 

9.14 

30.819 

141.69 

6 

CHERAPUNJI 

A1 

rs 

337.2 

L 

0.052 

21.47 

3.57 

15,839 

92.27 



A2 

33 

336.9 

T 

0.055 

26.78 

6.22 

15,219 

111 .90 



A3 

= 

277.0 

V 

0.024 

20.45 

3.81 

17,780 

65.50 

7 

DAUKI 

A1 

33 

308.5 

L 

0. 1 10 

46.87 

7,56 

15.779 

152.82 



A2 

33 

308.4 

T 

0.073 

44.53 

5.78 i 

14.899 

150.53 



A3 

== 

249.5 

V 

0.031 

20.74 

6.52 

18, 139 

73.07 

8 

DIPHU 

A1 

S 

190.1 

1 L 

0.283 

181.48 

23.25 

19,979 

636 , 62 



A2 

ss 

187.5 

1 T 

0.338 

205.55 

22.74 

16.779 

561 . 12 



A3 

= 

121.4 

1 ^ 

0.180 

56. 12 

9,00 

34.659 

196.14 

9 

DULOO 

A1 

zz 

232.6 

1 

1 

; L 

0.064 

) 

58.89 

14,40 

18.739 

274.51 



A2 

33 

233. 1 

T 

0-062 

53.49 : 

12.61 

20.099 i 

238 , 49 



A3 

- 

180.0 

V 

1 

0.039 

33.83 

7.63 

28.079 ! 

156,03 

10 

GUNJUNG 

A1 


211.2 

1 

1 

1 L 

0.092 

44. 19 

8.74 

28.199 ! 

176.73 



A2 

33 

210.7 

1 T 

0. 133 : 

52,05 

8.58 

25,239 

202 . 44 



A3 

= 

150.4 

1 

0.062 

25.56 

6.74 

41-279 

88.27 

1 1 

HAJADISA 

A1 

= 

182.5 

1 

1 ^ 

0.092 

42.73 

9.72 

34,180 

175. 12 



A2 

33 

181.6 

T 

0.098 

46.43 

8.77 

32.180 

174,97 



A3 

= 

120.1 

V 

0.046 

22.28 

5.98 

43.599 

89, 18 

I 12 1 

HARENGJAO 

A1 

33 

225.0 

L 

0.065 

40.99 

7.57 

20.839 

167.89 

1 

1 


A2 

33 

225.1 

T 

0.078 

44.39 

8.92 

19.459 

161. 06 



A3 

IS 

168.4 

V 

0.032 

22.15 

6.31 

21 .260 

69.89 


, . con td . 
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SR. 

NO. 

LOCATION 

EPl CENTRAL 
DISTANCE 

(Km) 

DIRN* 

PEAK 

6R0UND 

ACCL. 

(Q) 

PEAK 

GROUND 

VEL. 

<mm/s ) 

PEAK 

GROUND 

DISPL, 

(mm> 

DURATION 

TRIFUNAC 

BRADY 

<1975) 

(sec ) 

SPECTRUM 
INTENSITY 
FOR 5% 
DAMPING 

(mm ) 

13 

HOJAI 

A1 

= 

245.8 

L 

0.108 

51.88 

14,06 

22.439 

203.24 



A2 

r 

243.6 

T 

0. 134 

65.61 

11.61 

21.939 

260.09 



A3 

= 

177.5 

V 

0.061 

25.25 

7.83 

41.900 

121.04 

14 

JALLALPUR 

A1 

= 

265.2 

L 

0.029 

30,40 

4.81 

13.559 

130.38 



A2 

= 

265.5 

T 

0.023 

20.46 

5.18 

13.879 

82.49 



A3 

s 

209,8 

V 

0.015 

17.07 

4.01 

13.419 

57.62 

15 

JHIRIGHAT 

A1 

s: 

202,8 

L 

0,098 

88.30 

11.15 

18.500 

364.48 



A2 

= 

203-8 

T 

0.089 

70.35 

10.05 

20.779 

296.53 



A3 


155.3 

Y 

0.031 

30.38 

8.21 

25.700 

148. 15 

16 

KALAIN 

A1 

= 

253.8 

L 

0,057 

77.76 

15.87 

17.979 

328.96 



A2 

s 

254.2 

T 

0,051 

51.09 

15.71 

19.799 

223.70 



A3 

s; 

199.3 

V 

0.029 

43.21 

10.89 

22.620 

147.82 

17 

katakhal 

A1 

S5 

248,7 

L 

0.067 

110.74 

23,34 

18.360 

506.98 



A2 

S 

249,5 

T 

0,059 

1 93.56 

25,26 

19.339 

; 366.33 



A3 

= 

197.8 

V 

0.018 

1 33,72 

1 

10.11 

21 . 199 

1 133,94 

1 

18 

KHLIEHRIAT 

A1 


274.8 

L 

0,070 

' 29.21 

4,55 

26,780 

105.75 



A2 

= 

274.2 

T 

0,071 

31.74 

5.42 

31.659 

112.44 



A3 


213.4 

1 ^ 

1 

0,035 

1 

16.11 

3.44 

46.219 

1 47.66 

1 

19 

KOONBER 

A1 

— 

210.7 

1 

1 L 

0.049 

39.04 

9.86 

1 6 , 899 

184.96 



A2 

= 

211-2 

i T 

0.037 

34.88 

7.14 

18.920 

169.30 



A3 

= 

158.2 

V 

0.027 

I 

23.23 

1 

5.57 

19,920 

98.26 

20 

LOHARGHAT 

A1 

ss 

376. 1 

1 L 

0.058 

36,58 

8,58 

30,879 

205 . 47 



A2 

=: 

374.1 

i T 

1 0.054 

1 48,61 

9.91 

28,459 

255,16 



A3 


303.4 

I V 

0.022 

! 20,42 

5.06 

27,360 

j 

94.73 

21 

IMAWKYRWAT 

A1 

= 

364,6 

L 

0,046 

25,18 

3,10 

17,799 ! 

70,94 


1 

A2 

sz 

364, 1 

T 

0.046 

1 21,98 

4.80 

18,920 

78,82 



A3 

= 

303.4 

V 

0.032 

17,51 

4,69 

18.360 

66.87 

22 

NAWPHLANG 

A1 

s 

335.6 

L 

0.115 

‘ 39,17 

6.78 

25 . 359 

102.96 



A2 

= 

335.0 

T 

0.107 

38,99 

5.34 

29,659 

83,51 



A3 


273.4 

V 

0,036 

1 12.52{ 

! 

3,85 

30.019 

45.56 


contd 
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SR. 

NO. 

LOCATION 

EPICENTRAL 

DISTANCE 

(Km) 

DIRN* 

PEAK 

GROUND 

ACCL. 

(g) 

PEAK 

GROUND 

VEL. 

(mm/s) 

PEAK 
GROUND 
DISPL . 

<mm ) 

DURATION 
TRIFUNAC 
^ BRADY 
(1975) 

(sec ) 

SPECTRUM 
INTENSITY 
FOR 57. 
DAMPING 

(mm > 

23 

NAWSYNRAN 

A1 

= 

352.3 

L 

0.085 

33.67 

8,29 

17.399 

136.06 



A2 

= 

352.0 

T 

0.064 

28.61 

5.47 

18. 199 

124.34 



A3 

= 

292.0 

V 

0.036 

21.45 

5.89 

17,760 

90.86 

24 

NONGKHLAW 

A1 


352.8 

L 

0. 138 

64.50 

6.02 

30.879 

149,06 



A2 

= 

351 .8 

T 

0,146 

55.66 

11.07 

28.100 

182.58 



A3 

= 

280.4 

V 

0.083 

34.17 

9.46 

31 .399 

97.47 

25 

NONGSTOIN 

A1 


385.9 

L 

0.054 

24.09 

9.54 

28.399 ’ 

83,96 



A2 

= 

385.3 

T 

0.052 

31.04 

11.16 

27.779 

99.32 



A3 

= 

323.5 

V 

0.042 

14.80 

6.14 

29,540 

65.43 

26 

PANIHUR 

A1 


238.1 

L 

0.169 

55.19 

5.29 

31 .379 

132,72 



A2 

= 

236,7 

T 

0. 125 

44.38 

6.06 

41 .819 

110,47 



A3 

= 

172.3 

V 

0.073 

20.85 

6.14 

48.159 

59.68 

27 

PYNURSLA 

A1 

zz 

320,8 

L 

0-049 

42.13 

5.92 

22,439 

138.09 



A2 

= 

320.4 

T 

0.051 

32.24 

7.22 

24,200 

112,20 



A3 

s: 

260.2 

V 

0.036 

31.55 

i 

6.35 

23.159 

116,57 

28 

SA I ISAM A 

A1 


280.3 

L 

0.211 

87.31 

7.83 

33,799 

218.48 



A2 

= 

279,0 

T 

0.233 

103.63 

13.89 

32.519 

261,50 



A3 

= 

214.6 

V 

0.099 

33.77 

4.82 

46,840 

95.51 

29 

SHILLONB 

A1 

— 

323.9 

L 

0,075 

20.44 ! 

3.24 

24,319 

65.13 



A2 

=: 

323 . 1 

T 

0,057 

17,55 ; 

5.86 i 

26,139 

58.26 



A3 


260.5 

V 

0.036 

9.75 ! 

3.31 1 

26,000 

42.02 

30 

SILCHAR 

A1 

- 

233.0 

L 

0.064 

71,52 

21 , 10 

22,579 

277.62 



A2 

r; 

233.0 

T 

0.091 

98,76 

21 . 13 

18,079 

426,91 



A3 

= 

188. 1 

V 

0.025 i 

i 

29.98 

8,62 

41 ,720 

98.40 

31 

UMMULONG 

A1 


298.4 

i 

1 

0, 100 ! 

30.76 

7,86 

30,639 

87,60 



A2 

=: 

297.6 

T j 

0. 153i 

43,95 

12.80 

25.780 

128.93 



A3 


235.3 

V i 

0.062 i 

27.33 

11.28 

34.859 

88. 15 

32 

UMRONGSO 

A1 

= 

252.3 

L 

0.780 ! 

51.66 

6 . 46 

31,359 

197,65 



A2 

= 

251.4 

T 

0.081 

39.98 

7.76 

31 , 099 

169.49 



A3 


188.7 

V 

0.045 

17.48 

4.46 

34, 159 

51,01 1 

33 

UMSNING 

A1 

=: 

329.4 

L 

0.156 

41 ,46 

7.88 

29,159 

121 . 41 



A2 

=: 

328.2 

T 

0.153 

51 .94 

10. 12 

26,439 

149.39 



1 _ 

= 

264.3 

V 

0.071 

38.31 

12,50 

41 . 060 

61.92 


, T : Transverse, V ; Vertical 


■» L : Longitudinal 
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Table A -6 For October 20, 1991 earthquake recorded in UttarPradesh Array 


SR. 

NO. 

LOCATION 

EP I CENTRAL 
DISTANCE 

(km ) 

DIRNt 

PEAK 

GROUND 

ACCL. 

( 9 ) 

PEAK 

GROUND 

VEL. 

<mm/sec ) 

PEAK 

GROUND 

DISPL. 

(mm ) 

DURATION 
TRIFUNAC 
&( BRADY 

(1975) 
(sec ) 

SPECTRUM 
INTENSITY 
FOR 57* 
DAMPING 

(mm > 

1 

BH AT WAR I 

25.3 

L 

0.253 

178.73 

37.54 

13.559 

128.28 




T 

0.247 

297.78 

53.23 

11 .699 

206.27 




V 

0.294 

133.65 

23.53 

17.879 

77.09 

2 

UTTARKASHI 

40.2 

L 

0.242 

169.56 

21.15 

15.119 

207.06 




T 

0.310 

194.68 

19.85 

14.079 

172. 11 




V 

0.196 

141 .56 

22.98 

18.560 

134. 67 

3 

GHANSALI 

41.1 

L 

0.118 

80.44 

13.57 

11.119 

92.52 




T 

0.117 

78.21 

13.37 

11.139 

91.21 




V 

0.101 

95.94 

25.90 

13.079 

83.72 

4 

RUDRAPRAYAG 

54.4 

L 

0.053 

20.70 

7.85 

6.860 

506.97 




T 

0,052 

27.06 

4.01 

6.860 

579.37 




V 

0.045 

17.92 

3.83 

11.56 

433.27 

5 

TEHRI 

54.4 

L 

0 . 073 

42.15 

8.17 

7.820 

709.71 




T 

0.062 

92 . 30 

19.84 

1 5 . 460 

1 1382.31 




V 

0.059 

88.41 

123.68 

1 

7 . 660 

568.51 

6 

SRINASAR 

59.4 

L 

0.067 

19.45 

1 

' 5 . 86 

16.279 

361.40 




T 

0.050 

20.20 

5.07 

17.719 

345,22 




V 

0.034 

35.25 

7.62 

17.779 

395.03 

7 

BARKOT ^ 

62.7 

L 

0 . 095 

57.87 

10.93 

14. 159 

217.77 




T 

0 . 082 

44.84 

6.98 i 

13.539 

428.18 




V 

0.044 

27.53 

5.62 : 

14.299 : 

422.79 

8 

KOTESHWAR 

64.6 

L 

0.101 1 

51.55 

11.14 : 

14.140 

227.38 




T 

0.066 

39.27 

6.78 

1 8 . 330 

200.27 




V 

0.076 

85.25 

20.47 

19.520 

343.09 

9 

KARNAPRAYAG ’ 

65 . 5 

L 

0.062 

36.90 

5.80 

25.379 

111.96 




T 

0-079 i 

37.30 

4.02 

25.520 

93.24 




V 

0.026 ' 

14.98 

2.17 

27.719 

164.18 

10 

PUROLA 

76.0 

L 

0.075 

48.13 

8.44 

23.860 

77.60 




T 

0.093 

45.91 

9.22 

21 .739 

75,25 




V 

0.053 

25.57 

4.49 

22.639 

67-10 


- contd 
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SR- 

NO. 

LOCATION 

EPICENTRAL 

DISTANCE 

(km) 

DlRNt 

PEAK 
GROUND 
ACCL . 

( 9 ) 

PEAK 

GROUND 

VEL. 

1 ( mm/sec ) 

PEAK 
GROUND 
DIBPL . 

(mm ) 

DURATION 
TRIFUNAC 
h BRADY 

(1975 ) 
(sec ) 

SPECTRUM 
INTENSITV 
FOR SX 
DAMPING 

( mm > 

1 1 

jKOTI 

105.4 

L 

0.021 

23.43 

4.27 

12.920 

199.10 




T 

0.042 

28.60 

3.40 

1 14.019 

104.79 




V 

0.015 

1'7.65 

5-05 

16.859 

49.64 

12 

KOSANI 

143.6 

L 

0.029 I 

18.82 

3.77 

11.199 

54.75 




T 

0 . 032 ‘ 

15.55 

2.86 

1 1 . 439 

52.71 




V 

0.011 

9.17 i 

2.41 

12.139 

34.64 

13 

ALMORA 

150.5 

L 

0.018 

13.32 

3.42 

16.299 1 

45.40 




T 

0.021 

12.62 

4-50 : 

17.440 

46.78 




V 

0.019 

15.48 

3.98 

18.039 

44.25 


* L : Longitudinal, T : Transverse, V : Vertical 





